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Qh| Using techniques from loop quantum gravity, the standard theory of cosmological 

' perturbations was recently generalized to encompass the Planck era. We now apply 

OO ■ this framework to explore pre- inflationary dynamics. The framework enables us to 

isolate and resolve the true trans-Planckian difficulties, with interesting lessons both 
for theory and observations. Specifically, for a large class of initial conditions at the 
bounce, we are led to a self consistent extension of the inflationary paradigm over 
bi)' the 11 orders of magnitude in density and curvature, from the big bounce to the 

onset of slow roll. In addition, for a narrow window of initial conditions, there are 
. departures from the standard paradigm, with novel effects — such as a modification 

, of the consistency relation between the ratio of the tensor to scalar power spectrum 

I and the tensor spectral index, as well as a new source for non-Gaussianities — which 

could extend the reach of cosmological observations to the deep Planck regime of the 
early universe. 
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' I. INTRODUCTION 

>< 

I Thanks to a powerful interplay between theory and observations, our understanding of 

the early universe has advanced very significantly in recent years. Consequently, cosmology 
offers a natural arena to develop and test approaches to quantum gravity. 

In this paper we will focus on inflation, the leading scenario to successfully account for the 
tiny inhomogeneities observed in the cosmic microwave background (CMB) (see, e.g., [1-5]). 
An impressive feature of this scenario is that it involves a rather small set of assumptions. 
As noted in [6, 7], these can be summarized as follows: 

• Sometime in its early history, the universe underwent a phase of rapid expansion during 
which the Hubble parameter was nearly constant; 

• During this phase, the universe was well described by a Friedman, Lemaitre, Robertson 



*Electronic address: I.AgulloRodenas@damtp.cain.ac.uk 
^Electronic address: aslitekar@gravity.psu.edu 
•'•Electronic address: nelson@gravity.psu.edu 



2 



Walker (FLRW) solution to Einstein's equations with a scalar field as matter source, 
together with small inhomogeneities which are well approximated by first order per- 
turbations; 

• Fourier modes of the quantum fields representing perturbations were initially in the 
Bunch-Davies (BD) vacuum at least for a certain set of wave numbers^; and, 

• Soon after any mode exits the Hubble radius, its quantum fluctuation can be regarded 
as a classical perturbation and evolved via linearized Einstein's equations. 

Analysis of these inflationary perturbations implies that there must be tiny inhomo- 
geneities at the surface of last scattering whose detailed features have been observed in the 
CMB. Furthermore, these tiny inhomogeneities serve as seeds that grow into a large scale 
structure which is in excellent qualitative agreement with observations. Therefore, even 
though the underlying assumptions are by no means 'obvious' or 'compelling', the success 
of the subsequent predictions is impressive. 

Yet, as is well known (see e.g. [8]), the scenario is conceptually incomplete in several 
respects. In particular, as Borde, Guth and Vilenkin [9] showed, inflationary space-times 
inherit the big-bang singularity. Physically, this occurs because one continues to use general 
relativity even in the Planck regime in which it is simply not applicable. It is widely expected 
that new physics in this regime will resolve the singularity, significantly changing the very 
early history of the universe. Will inflation arise naturally in the resulting deeper theory? 
Or, more modestly, can one at least obtain a consistent quantum gravity extension of this 
scenario? 

To fully address this question one has to face two distinct sets of issues. The first set 
has its origin in particle physics. What is the physical origin of the scalar field that plays 
the role of the inflaton? It should arise naturally in the fundamental unified theory. Can 
we derive the inflationary potential from some first principles? Is there a single inflaton, 
as in the simplest models that have been successful so far, or many? How do they interact 
with one another? Even more importantly, what are the interactions that produce the 
elementary constituents of the standard model of particle physics during reheating? These 
issues remain open, but fall in the realm of theories aimed at unification. In this and the 
companion papers, we do not address them. 

Rather, we will continue to use the simplest, single inflaton model and focus on the second 
set of issues related to gravity. In the systematic evolution starting from the Planck regime, 
does a slow roll phase compatible with the WMAP data arise generically in the background 
geometry, or, is an enormous fine tuning needed? In classical general relativity, if we evolve 
the modes of observational interest back in time, they become trans-Planckian at some stage. 
Is there a quantum field theory (QFT) on quantum space-times that is needed to adequately 
handle physics at that stage? Can one arrive at the BD vacuum at the onset of the slow 
roll from more fundamental considerations? Or, is an even more elaborate fine tuning of 
the quantum state of perturbations necessary in the Planck regime? The natural state 
of quantum perturbations at the onset of inflation may be indistinguishable from the BD 
vacuum for currently available data but, because of the pre- inflationary dynamics especially 



More precisely, if fco denotes the co-moving wave number of the mode which has just re-entered the Hubble 
radius now then, at the time t{ko) at which this mode exited the Hubble radius during inflation, Fourier 
modes with co- moving wave numbers in the range ^ (fco, 2000fco) were in the Bunch-Davies vacuum. 
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in the Planck regime, it could well carry certain excitations over the BD vacuum. These can 
then source non-Gaussianities for the subsequent evolution during inflation [10-14] which, 
in turn, could give rise to novel effects in the CMB and galaxy distribution [15-17] which 
could be observed in future missions. Can the state at the onset of the slow roll be close 
enough to the BD vacuum to agree with current observations and yet be sufficiently different 
to give rise to such effects? 

This is the third in a series of papers whose goal is to address these questions that origi- 
nate in quantum gravity. The first paper [18] presented a succinct summary of underlying 
framework and the principal observational consequences for a broad audience. The second 
[19] was primarily addressed to the quantum gravity community and focused on conceptual 
and mathematical aspects of the framework. In this paper we use that framework to an- 
alyze phenomenological issues of interest to cosmologists within the standard infiationary 
scenario. 

The main ideas can be summarized as follows. 

By now there is a large body of results in loop quantum cosmology (LQC) that provides 
a detailed description of the Planck scale physics in a number of cosmological models. These 
include the k=0 and k=l (FLRW) models [20-28], possibly with a non-zero cosmological 
constant [29-31], the anisotropic Bianchi I, II and IX models [32-35] and the simplest of the 
inhomogeneous models — the Gowdy space-times — widely studied in exact general relativity 
[36-40]. In all cases, the big bang singularity is resolved and replaced by quantum bounces 
(for a review, see, e.g., [7, 41]). Thus, the first conceptual limitation of the standard general 
relativistic models is overcome. It is therefore natural to use LQC as the point of departure 
for extending the cosmological perturbation theory. How do new effects in the Planck 
domain — such as the robust superinfiation phase that follows the density bounce — alter 
the standard infiationary scenario? Do generic initial conditions at the bounce lead to the 
desired slow roll infiation later, or, is a fine tuning necessary? 

To incorporate cosmological perturbations, however, one cannot just mimic the standard 
procedure used in general relativity: Since loop quantum gravity (LQG) is still incomplete, 
we do not yet have the analog of full Einstein's equations to perturb. A number of LQC 
inspired strategies have been devised to overcome this obstacle.^ In this series of papers 
we follow a mainstream strategy in LQG: First truncate the classical theory in a manner 
appropriate to the physical problem under consideration, then carry out quantization paying 
due attention to the underlying quantum geometry of LQG, and finally work out the con- 
sequences of the resulting framework. This procedure has been successful not only in the 
cosmological models referred to above but also in the investigation of quantum black holes 
in LQG [42-44] and derivation of the graviton propagator using spin foams [45, 46]. 

For the present application, we will truncate the full phase space, keeping only the FLRW 
backgrounds with first order inhomogeneous perturbations. To facilitate easy comparison 
with the standard quantum theory of cosmological perturbations, we will employ a hybrid 
scheme used in [36-40, 47], in which one uses the standard LQC for homogeneous modes 
and 'Fock-like' states for the first order perturbations. There is however a key difference 
from the standard cosmological perturbation theory: Because the background space-time 
now has quantum geometry, we are led to use QFT on quantum cosmological space-times. 
This theory was introduced in [48] emphasizing conceptual issues and developed in [19] to 
incorporate the infinitely many modes of perturbations and the associated regularization 
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and renormalization techniques which, as discussed below, are key to checking the self- 
consistency of our initial truncation. Since we only have a probability amplitude for various 
FLRW geometries — rather than a sharply defined classical FLRW metric — at a fundamental 
level evolution has to be described using a relational time variable rather than cosmic or 
conformal time. Thus, quantum gravity introduces genuinely new conceptual elements. But 
one can systematically 'descend' to a description in terms of cosmic and conformal time of 
an effective metric that incorporates the quantum gravity corrections. 

Since quantum perturbations now propagate on quantum geometries which are all regular, 
free of singularities, the framework automatically encompasses the Planck regime. What is 
the status of the 'trans-Planckian issue' which, in heuristic discussions, is often associated 
with modes of trans-Planckian frequencies? Recall, first, that the quantum Riemannian 
geometry underlying LQG is subtle [49, 50]. For example, while there is a minimum non- 
zero eigenvalue of the area operators, there is no such minimum for the volume operators, 
even though their eigenvalues are also discrete [51, 52].^ As a consequence, there is no 
fundamental obstacle preventing trans-Planckian modes of perturbations in our truncated 
theory. Indeed, in the homogeneous LQG models that have been analyzed in detail, the 
momentum p(^^) of the scalar field is generally huge in Planck units. This poses no problem 
and, in particular, on the physical Hilbert space the total energy density is still guaranteed to 
be bounded by 0.41ppi where ppi is the Planck density. Similarly, for perturbations, there is 
no a priori difficulty with trans-Planckian momenta (or frequencies) in the truncated theory 
considered here. The real danger is rather that, in presence of such modes, the energy density 
in perturbations may fail to be negligible compared to that in the quantum background 
geometry. This is a very non-trivial issue especially in the Planck regime following the 
bounce and requires a careful treatment of regularization and renormalization of the stress 
energy tensor of quantum perturbations. If the energy density does become comparable to 
that in the background, then we would not be able to neglect the back-reaction and our 
truncation would fail to be self-consistent. This is the trans-Planckian problem we face in 
our theory of quantum perturbations on inflationary quantum geometries. 

The paper is organized as follows. Section II summarizes the underlying conceptual and 
mathematical framework (constructed in [19] using prior results from [6, 22-25, 56, 57]) for 
the cosmology community. This discussion encompasses the scalar as well as the tensor 
modes. Section III discusses the initial conditions for the wave function \I^o representing the 
background quantum geometry, and for the wave function tp representing the perturbations 
propagating on the quantum geometry \l/o. The choices can be motivated by general sym- 
metry principles and physical considerations. Sections IV and V present results of numerical 
simulations of the pre-inflationary dynamics for the quantum corrected background and per- 
turbations. The free parameter that dictates the evolution of the background is the value 
0B of the inflaton at the bounce. In section V we vary this parameter (within a range that 
is numerically feasible) and explore consequences on the observed power spectrum at the 
end of inflation. We flnd an unforeseen and interesting structure. The LQC power spectrum 
is essentially independent of the value of 0b ■ But recall that there is only a flnite inter- 
val in the space of co-moving wave numbers k that is relevant for the CMB observations. 
This interval is completely determined by the value of (pB, moving steadily to the right along 
the k-axis as 0b increases. If 0b is sufficiently low, the quantum state at the onset of the 



^ Properties of the eigenvalues of length operators [53-55] have not been analyzed in comparable detail. 
But since their definitions involve volume operators, it is expected that there would be no 'length gap'. 
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slow roll is sufficiently different from the BD vacuum to lead to observable effects for future 
measurements along the lines of [15-17]. 

In section VI we analyze the issue of self consistency of the truncation strategy: Does 
the stress energy in the perturbations remain small compared to that in the background 
throughout the evolution from the big bounce to the onset of the slow roll? The main result 
is that the self-consistency criterion is met for a large class of initial conditions. In these 
cases, our main approximation — the initial truncation — is viable. Had it failed, it would 
have been inconsistent to keep only the first order perturbations. Then one would have 
had to wait for significant advances in full LQG to extend the inflationary paradigm to the 
quantum gravity regime. Thus, technically as well as conceptually it is quite non-trivial that 
this does not happen and one obtains a self consistent extension of the standard inflationary 
paradigm all the way to the deep Planck regime of the big bounce for a wide range of initial 
conditions. 

Of course, self-consistency by itself does not imply that our truncated solution is neces- 
sarily close to an exact one (because the sum of all higher order terms could be large). But 
this limitation is common to all perturbation theories, classical or quantum. In particular, in 
classical cosmology the total back reaction is routinely neglected if stress-energy in the ffist 
order perturbations is small compared to that in the background, even though strictly this 
does not imply that there is an exact solution close to the linearized one. If the first order 
perturbations do not grow so much in time as to become comparable to the background, 
the test field approximation is self consistent and is regarded as trustworthy. We adopt the 
same viewpoint here. 

In section VII we summarize the main results and discuss open issues. The Appendix 
summarizes the technical differences between Ref. [19] and the underlying framework pre- 
sented in section II which is tailored to the inflation. 

We will use the following conventions. The signature of the space-time metric will be 
-,-!-,+,+. We set c=l but keep G and h explicitly in various equations to facilitate the 
distinction between classical and quantum effects. As is usual, we will set k = SnG and, 
for the Barbero Immirzi parameter 7 of LQG, use the value 7 ~ 0.24 that comes from 
the black hole entropy calculations. Finally, we will use Planck units used in the quantum 
gravity literature rather than the reduced Planck units often used in cosmology. (Thus, our 
Planck mass mpi = \fh[G is related to the reduced Planck mass Mpi via mpi = a/SttMpi.) 
Numerical values of all quantities are given in dimensionless Planck units ipi = mpi = t-pi = 
1. 

II. THE UNDERLYING FRAMEWORK 

In discussions of the early universe, one generally begins with a FLRW solution to Ein- 
stein's equation with suitable matter fields as sources, then considers the space of solutions to 
the ffist order perturbation equations on this background space-time, and finally quantizes 
them using standard techniques of QFT in curved space-times. To incorporate the Planck 
regime, the FLRW background must also be treated quantum mechanically. Therefore, 
we cannot begin with a solution to Einstein's equations. Instead, we will ffist introduce 
a Hamiltonian framework that encompasses both the FLRW backgrounds and ffist order 
perturbations thereon and then pass to the quantum theory of the combined system as a 
whole. As explained in section I, this will lead us to a quantum theory of fields representing 
linearized perturbations, propagating on a quantum FLRW geometry. 
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We will assume that the spatial topology is and focus on the k=0 FLRW case. As 
in the standard inflationary models, the matter fleld will be taken to be a scalar fleld and 
for detailed calculations we will use the simplest potential V{(f)) = (l/2)m^0^. Most of the 
results discussed in this section are generalizations of those discussed in detail in [19], now 
allowing for a inflaton potential V{(j)) and M"^ spatial topology (in place of a 3-torus). 



A. The truncated phase space 

Let us flrst truncate the full phase space of general relativity to the physical problem 
of interest. This truncated phase space FTrun is of the form Fxrun _= To x Ti where To 
is the 4-dimensional phase space for the FLRW backgrounds and Fi of gauge invariant 
perturbations. 

Because the background flelds are homogeneous and the spatial topology is M.^, the spatial 
integrals that appear in the expressions of the symplectic structure and Hamiltonians on Fo 
diverge. Therefore, in the construction of the Hamiltonian framework of the homogeneous 
sector, one flrst flxes co-moving Cartesian coordinates x", introduces an elementary cell C 
whose edges are aligned along these coordinates and have equal coordinate lengths £, and 
restricts all integrations to C. This is an infrared cut-off for the homogenous sector which is 
removed at the end by letting i tend to inflnity [7]. 

In LQC, in place of the scale factor a and its conjugate momentum 7i(^a), it is customary 
to use the following pair: 



Aa^f and b - 



'^=— r; and b = -^^7r(,) (2.1) 



so that the fundamental Poisson bracket is given by {b, u} = j^. Thus, apart from constants, 
u gives the physical volume of the flducial cell C and, on solutions to Einstein's equations, 
b is related to the Hubble parameter H = a/a via b = 'yH. The homogeneous scalar fleld 
and its momentum will be denoted by and respectively. Thus the background FLRW 
phase space Fq is coordinatized by (z^, b; and carries a single scalar (or Hamiltonian) 

constraint: 



So[A'hom] — A'hom 



47 K'-fhu 4 



0, (2.2) 



where Ahom represents a homogeneous lapse. Einstein dynamics is generated by this scalar 
constraint. The choice A^hom = 1 yields evolution in cosmic time; Ahom = a in conformal 
time; and Ahom = cl^ in harmonic time. (A^hom = 0?^^ /P{4>) corresponds to using the inflaton 
(f) as time in 'patches' of dynamical trajectories along which it is monotonic.) 

The phase space Fi of gauge invariant perturbations is spanned by 3 canonically conjugate 
pairs, one representing the Mukhanov- Sasaki scalar mode and two representing the tensor 
modes. It is simplest to work in the (co-moving) momentum space and represent them by 

^Qk^'^f^'^P-^ pf'P?'P?)- (2-3) 

f\i t\i r\j r\j IX 

Since the perturbations are not homogeneous, in the discussion of Fi one can work directly 
with the entire M^; restriction to the cell C is not necessary. We choose this route to 
avoid an artiflcial quantization of the momenta k that the restriction to the cell would have 
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introduced. Technically, we assume that the perturbations are square-integrable so that one 
can freely pass between x and k spaces. Then the Poisson bracket relations are given by 
{Qk^ P^p''} ~ (2vr)^(5(A; + k'), and similarly for tensor modes. For simplicity of notation, we 

will often use collective labels 7^, for the two tensor modes. 

On the truncated phase space Fxrun, the dynamical trajectories follow integral curves of 
the 'evolution vector field' ^Dyn = ^o^djjEio + f2"^(9/3§2 where flo is the symplectic structure 
on Fq, Qi on Fi and So is defined in Eq (2.2). The function §2 governing the time evolution of 
perturbations is obtained from the second order truncation §2 of the full scalar constraint of 
general relativity by keeping only the terms which are quadratic in first order perturbations 
and deleting the term that is linear in the second order perturbations.^ It is given by 



and 




In (2.5), 21 is a function only of the background variables (and therefore has no k- 
dependence): 

01= {fV {(!))- + (2.6) 

with f = {3Kpf^^) /{{l/2)pf^^+i^a^V{(j))) (see Appendix (A). Since p(^) = a^i^<j), the quantity 
f is essentially the kinetic fraction of the total energy density.) 

The resulting time evolution of tensor modes is the same as that of a massless scalar field 
on a given FLRW background. By contrast, the scalar mode experiences an external, time 
dependent potential 21 in addition to the FLRW geometry. If V{(j)) were to vanish, we would 
have 21 = and then the dynamics of the (trivially rescaled) scalar mode Q^/2a/k would be 
the same as that of the tensor mode 7^ (as in [19]). 

Recall first that the time evolution in the FLRW sector is generated by the restriction So 
to the FLRW subspace Fq of the Hamiltonian constraint on the full phase space F of general 
relativity. This evolution is reproduced by the first term Q'^^di^Eio of the dynamical vector 
field ^Dyn- Consequently these dynamical trajectories are solutions to the full Einstein 
equations. However in the second part of the dynamical vector field, Sj is not constrained 
to vanish. Indeed, even in the second order truncation of the full constraint on F, it is S2 
that vanishes, not Sg. The fact that dynamics is governed by a constraint only on Fq and 
not on the full Fxrun has an important consequence in the quantum theory. 



^ Here the Greek indices refer to the (infinite dimensional) tangent space to FTrun- In terms of Poisson 
brackets, the evolution of functions F of background quantities z^, b, is given hy F = {F, So}o) 

where the 'dot' denotes derivative with respect to the time variable determined by the lapse and the 
Poisson bracket is taken only with respect to the background variables. The evolution of functions / of 
gauge invariant perturbations (Qt,Tt', P-^\p-^^) is given by / = {/, Snii where the Poisson bracket is 
taken only with respect to the perturbations, treating the background quantities in the expression of §2 
(and possibly /) as given time- dependent variables determined by first solving the evolution equation for 
functions F. 
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Finally, because §2 depends not only on perturbations but also on the background variable 
V (i.e. the scale factor), the dynamical vector field is not of the form i^l^ + fi"'^) d^H for 
any function H on FTmn: In striking contrast to full general relativity, or its FLRW sector, 
dynamics is not generated by any Hamiltonian on Fttuii- In the phase space language, we 
have to first obtain the dynamical trajectory on Fq and then lijt it to the truncated phase 
space Fxrun using integral curves of X^^^. This procedure just refiects the steps normally 
followed using space of solutions: one first solves for the background scale factor 0(77) and 
the scalar field 0(77) using FLRW equations, fixes a specific solution as the background 
space-time and then solves for perturbation equations on this background. Therefore, in the 
quantum theory we are led to first solve for the background wave function \l/o(^, 4>) and then 
'lift this quantum trajectory' to a wave function \l/o (g) i'iQj:, T^, 4>) describing the evolution 
of the quantum wave function on the quantum background geometry \l/o. 

B. Quantum FLRW geometry in LQC 

Since the phase space of the truncated system is a product, Fxrun = To x Ti, the Hilbert 
space of quantum states has the form "H = 'Ho^'Hi- In this sub-section we will focus on the 
space Tio of quantum states of the background geometry and the metric operator thereon. 
The sector "Hi of perturbations will be discussed in the next sub-section. 

1. Quantum dynamics of the FLRW background 

Although we have used background variables tailored to LQC, the classical Hamiltonian 
theory is the same as in general relativity. It is in the passage to the quantum theory that we 
use the LQG techniques. For simplicity, we will first summarize the basic ideas and return 
to two important technical points at the end. 

Recall that Fq is spanned by two canonically conjugate pairs (z/, b,p((^)), and carries a 
scalar constraint Sol^^^b) = (see (2.2)). Therefore quantum states are represented by wave 
functions ^o(i^5 4>) and the Dirac quantization procedure would lead us to impose So^o = 0, 
which takes the form {h'^d'^ + 0) \I^o = for a specific operator G. A careful analysis [58-60] 
of constrained systems implies that the physical Hilbert space Tio is spanned by states that 
satisfy a mathematically 'sharper' Hamiltonian constraint, which can be intuitively thought 
of as the 'positive frequency square-root' of Sq^o = [23]: 

-ihd^^oii^, 0) = Ho^oii^, 0) (2.7) 

where Hq is a self-adjoint operator whose explicit form will not be needed in this summary. 
Thus, the constraint is 'de-parameterized': its form suggests that we can interpret the scalar 
field as an 'internal' or a 'relational' time variable (a la Leibniz) with respect to which the 
'true' dynamical variable z> (and also the perturbations Q^, 7^) evolve. Thus, in the physical 
sector of the theory, is just a parameter; there is no longer an operator associated with it. 
Now, given \l/o(i^; 0b) at the bounce time = 0b, Eq (2.7) enables one to 'evolve' this initial 
state to obtain the quantum state "^oi^^, 0) at all times 0. This quantum evolution is non- 
singular; the wave function undergoes a bounce when the density reaches a maximum value 
Pmax ~ 0.41ppi. How does LQC evade the standard singularity theorems? The singularity 
theorem due to Borde, Guth and Vilenkin [9], tailored to infiation, is evaded because the 
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LQC universe has a contacting phase in the past, violating their assumption of eternal 
expansion. And LQC bypasses the original singularity theorems [61] in general relativity 
even in cases when matter satisfies all energy conditions because quantum geometry effects 
modify the (geometrical) left hand side of Einstein's equations. 

In the classical theory, a FLRW solution corresponds to a trajectory (f){t), or, elimi- 
nating the parameter t, just z/(0). A solution "^oii^, 0) to (2.7) represents the quantum analog 
of this trajectory. Since there is no sharp trajectory, we no longer have a single, sharply 
defined space-time metric in the background. Consequently, we do not have a canonical 
parameter representing the cosmic or conformal time. At a fundamental level, all dynamics 
in the Planck regime refers to the relational time (p. 

We will be interested in a specific class of states \E'o- To specify that class, it is convenient 
to first set V{(j)) = to get a qualitative insight into LQC dynamics. This case has been 
analyzed in great detail in the LQC literature. One can now start at late times, when the 
classical approximation is clearly excellent, and consider wave functions which are sharply 
peaked at a point on the phase space. Under evolution, these states are known to remain 
sharply peaked and the peaks are known to follow trajectories satisfying certain effective 
equations that encode the leading-order quantum corrections. These effective solutions are 
in excellent agreement with the FLRW trajectories of general relativity when the matter 
density and the curvature are less than a thousandth of the Planck scale [7, 62, 63]. But 
they strongly deviate from general relativity in the Planck regime: they define singularity- 
free, bouncing trajectories. This behavior has been established analytically in the k=0, 
A=0 case [60] and was checked numerically in the k=l, as well as A ^0 cases [27, 29, 31]. 
The effective equations themselves are rather general and continue to be meaningful when 
V{(p) 7^ 0, or if we replace the scalar field by other matter sources [64, 65]. 

In presence of the quadratic potential V{(j)) = (l/2)m^<^^, sharply peaked wave functions 
have been constructed whenever numerical simulations are feasible, which corresponds to 
the case when the kinetic energy is dominant at the bounce [56].^ Again, the peaks of these 
wave functions follow solutions to effective equations. We do not see any a priori reason why 
this general behavior will not continue away from kinetic energy dominated bounces. More 
importantly, as we will see in sections IV and V, it turns out that states \E'o which undergo a 
kinetic energy dominated bounce are the most interesting ones for potentially new physics. 
Therefore, while our underlying framework is valid for all states, in our detailed analysis and 
numerical simulations we will restrict ourselves to those solutions "^oi^^, 4>) to (2.7) which are 
sharply peaked on an effective trajectory from the bounce (at least) until the trajectory enters 
the regime in which the matter density and curvature are so low that general relativity is an 
excellent approximation. 

Since the FLRW sector of the phase space is a system with a Hamiltonian constraint, 
we used the Dirac quantization procedure and reinterpreted the constraint as an evolution 
equation with respect to the internal time 0. This is why we were naturally led to the 
Schrodinger picture in which states evolve and operators don't. However, in the cosmology 
literature one generally uses the Heisenberg picture. The transition can be carried out in 
the usual fashion. In the Heisenberg picture, states are frozen in time and it is the volume 
operator (or scale factor operator a = {n'jh |z>|/4£^)-^/^) that evolves with respect to (j) and 



^ In this case, it takes about 10^ Planck seconds to reach the onset of slow roll starting from the bounce 
(see, e.g., Table 1). As the potential energy at the bounce increases, this period increases rapidly and it 
has been difficult to maintain the required accuracy for these longer periods. 
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the space-time metric operator is given by [19]: 

gabdx^dx^ = = H;^ f d'^{(f)) H;U(j)'^ + d\(f)) dx^ (2.8) 

where we have used the fact that, in the classical theory, the lapse corresponding to the 
scalar field time is iVhom = 0'^^^/P{,f,), and, on physical quantum states, = —ihd^ = Ho- 
The operator Qab lives on a 4-manifold M coordinatized by 0, x, where = const surfaces are 
regarded as 'space' sections. The effective trajectories correspond to using the expectation 
value a = {^od'^o) in place of d in (2.8). 

At a fundamental level, then, the parameter serves as the time variable and the 
quantum geometry is determined by the physical (Heisenberg) operators a(0). However, as 
we will see in the next sub-section, an unforeseen simplification occurs in the dynamics of 
perturbations even in the Planck regime, enabling us to cast this exact quantum dynamics 
in terms of the conformal time of a quantum corrected metric. 

Remark: Mathematically as well as physically, the deparametrization procedure outlined 
above is straightforward in the case when the scalar field potential V{(f)) vanishes (for details, 
see [7, 23, 24]). In the present case one has to incorporate some subtleties on both fronts 
[56]. 

On the mathematical side, it turns out that the formal quantum constraint So again has 
the operator form h'^d'^ + O as in the V{(j)) = case but now the operator 6 is no longer 
essentially self-adjoint [7]. Therefore, to arrive at (2.7), one has to choose a self-adjoint 
extension of 0, and Ho is then the square-root of the corresponding operator |G)| which is 
self-adjoint by construction. This overall situation is the same as in the case when V{(j)) = 
but there is a positive cosmological constant A. That case is well-understood [31] and, for 
states of physical interest, the results have been shown to be quite insensitive of the choice of 
the self-adjoint extension. In the present case, the existing results indicate that the situation 
should be similar. 

On the physical side, in the V{(f)) = case serves as a global time parameter along 
classical dynamical trajectories. This is no longer the case when V{(j)) ^ 0. Consequently, 
one has to work with 'patches' of dynamical trajectories such that is monotonic in each 
patch. In simple examples, the associated subtleties in the quantum theory have been 
discussed in the literature (see, e.g., [66, 67]). But for our purposes, the situation is simpler 
because we can work just in one 'patch'. This is because we will start at the bounce with 
increasing and by the time it 'turns around' after climbing up the potential, the energy 
density and curvature are so low compared to the Planck scale that general relativity is a 
good approximation. Thus, is in fact single valued in the regime in which full quantum 
treatment of the background geometry is needed. 

2. Effective equations 

The effective equations of LQC, which track the evolution of the peaks of wave functions 
\E'o, have a number of consequences that seem surprising if one's intuition is based largely 
on inflationary dynamics within general relativity. In this sub-section we will summarize 
these features. While our focus is on the quadratic potential, most of these results hold for 
general (regular) potentials. 
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• The Hubble parameter: Recall that on solutions to Einstein's equations, the momentum 
b conjugate to u is related to the standard Hubble parameter H = a/a via b = '-/H. However, 
on solutions to the LQC effective equations we have 

H = ^ sin2Ab ^ (0.93) sin2Ab, (2.9) 
27A 

where A^ = Ay/Sivjipi ~ 5.2 is the 'area-gap' that sets the discreteness scale of LQC. Thus, 
in striking contrast to general relativity, the Hubble parameter H is bounded in LQC by 0.93 
in Planck units, b ranges over (0, vr/A); where b = vr/A corresponds to the infinite past and 
b = to the infinite future. General relativity is recovered in the limit A — ?■ 0, i.e., when 
quantum geometry effects can be ignored. In this limit, h ranges over (0, 00) and b = 00 
corresponds to the big-bang. Finally, note that while H is monotonic in general relativity 
(in the absence of spatial curvature and provided p > 0), (2.9) implies that this is not the 
case in LQC. 

• The Friedmann equation: Quantum geometry corrections modify the left side of Ein- 
stein's equations. In particular, the Friedmann equation becomes 

sin^Ab SnG SttG , , 

= p = (— + VU)). 2.10 

To compare with the standard Friedmann equation = (SnG/S) p, it is often convenient 
to write (2.10) using (2.9): 

where Pmax = S/SttG'j'^X^ ~ 0.41ppi. Thus the familiar linear relation between the Hubble 
parameter and the matter density no longer holds. Nonetheless, as is obvious from Eqs (2.9) 
- (2.11), away from the Planck regime — i.e., when A6 <^ 1, or, p <^ Pmax — we recover 
classical general relativity. 

• Bounces: In general relativity, the Friedmann equation implies that, if the matter 
density is positive, d cannot vanish (in the absence of curvature or a cosmological constant). 
Consequently, in any given solution, the universe is either contracting, or an expanding. 
By contrast, (2.11) now implies that z> vanishes at p = Pmax! this is the quantum bounce. 
To its past, the solution represents a contracting universe with z> < and to its future, an 
expanding one with z> > 0. 

In LQC, b is monotonically non-increasing (if we ignore the exceptional de Sitter solutions 
that exist for certain potentials), evolving from b = vr/A in the infinite past to in the 
infinite future. Eqs (2.10) and (2.11) imply that b = 7r/2A at the bounce. Thus, each 
solution undergoes precisely one bounce. 

• Boundedness of physical quantities: Intuitively one can think of the big bounce of LQC 
as the replacement of the big bang in general relativity. However, all physical quantities 
remain bounded at the big bounce. 

As we have already noted, the matter density achieves its maximum value pmax at the 
bounce. When the potential V{(f)) is bounded below, the Ricci scalar — the only non-trivial 
curvature scalar in these models — is bounded above by ~ 62 and is bounded above by 
^ 10.29. If the potential V{(j)) is bounded below, say V > Vo, then it follows from (2.10) 
that 0^ is bounded by 2pmax — 2K- If V grows unboundedly for large |0|, then 101 is also 
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bounded. For example, for V = m^0^/2, we have m|0|jnax ~ 0.90. Finally, one can also 
show that if z/ 7^ initially, it cannot vanish in finite proper time along any solution. 

Perhaps the most striking contrast with general relativity occurs in the behavior of the 
Hubble parameter H: It vanishes at the bounce while in general relativity it diverges at the 
singularity and is large throughout the Planck regime. As one would expect, H is bounded 
above in LQC, \H\ < 0.93, and achieves its upper bound in every solution at the end of 
super-inflation (b = tt/AX) that follows immediately after the bounce. 

To summarize, solutions to the effective equations in LQC are everywhere regular irre- 
spective of whether one focuses on matter density, curvature or the scale factor. 

C. QFT on quantum FLRW space-times 

In striking contrast to the FLRW solutions to Einstein's equations, the quantum geometry 
\E'o is regular; there is no big bang singularity. The quantum perturbations Q^, 7^ propagate 
on this regular quantum geometry. Therefore this framework is well suited to incorporate 
Planck scale physics by facing various quantum gravity issues directly. Our task is to specify 
the dynamics of quantum fields Qj:, 7^ on any given quantum background \&o(^, 0)- At first 
this task seems formidable. However, the detailed framework of QFT on quantum FLRW 
space-times introduced in [48] and further developed in [19] shows that there is an unforeseen 
simplification. 

1. Dressed, effective metric 

Let us start with the tensor perturbations. In the classical theory, the dynamics of 7^ is 
the same as that of a zero rest mass scalar field on the background FLRW metric. In the 
quantum theory, we have the following surprising result. 

• Within the test field approximation inherent to our truncation strategy, the dynamics 
of Tj: on any one of our background quantum geometries is completely equivalent to 
that of the quantum field 7^ propagating on a smooth but quantum corrected metric 
gab given by 

gabdx°'dx^ = ds^ 

where 

Here the expectation value is taken in the quantum geometry state and a(0) is the 
(Heisenberg) operator on Ho describing the scale factor. This is an exact mathematical 
equivalence. 

Let us now turn to the scalar perturbations. In the classical theory, the dynamics of 
is governed not only by the background FLRW metric g^h but also by a time dependent 
external potential 21 (see 2.5). In the quantum theory we have: 



r a^(0) + dx' . 



(2.12) 



and 



d' 



{Ho~'d\(t))Ho~') 



(2.13) 
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• Again, within the test field approximation, the dynamics of Qj: on any one of our 
background quantum geometries \l/o is completely equivalent to that of the quantum 
field Qj: propagating on a smooth but quantum corrected geometry gab given by (2.12) 
and a quantum corrected time-dependent potential 21, given by 

^^^^^(H;hHmi<P)aH<P)HP)_ (2.14) 

For tensor modes, this result was derived systematically in [19] and the generalization to 
the scalar modes follows the same steps. The essence of the result is that while the quantum 
perturbations T{x, (f) and Q{x, 0) propagate on the quantum geometry of they do not 
experience all the details of quantum fluctuations encoded in this state. The net effect of 
the quantum geometry on their dynamics^ is neatly encoded in just three time- dependent 
quantities constructed from ^o'- P{(f)), o,, 21. This is similar to what happens when light 
propagates in a medium. While it interacts with the atoms of the medium, the net effect 
can be encoded in just a few parameters such as the refractive index. Quantum geometry 
can be thought of as the medium through which the quantum fields Q, T propagate and the 
net effect of the medium is to provide quantum corrections (2.13) and (2.14) to t], a and 2t. 

Next, note that in this exact result the effective quantities that enter are not just ex- 
pectation values of the corresponding operators; they are much more sophisticated and, in 
particular, they are sensitive also the quantum fluctuations in \l/o. Indeed their definitions 
are quite intricate and could not have been arrived at without a detailed calculation. In 
particular, in the classical theory P((f,){(f>) and a^cj)) not only determine the metric gab but also 
the external, time dependent potential 21 experienced by the scalar perturbation Qg. In the 
quantum corrected description, by contrast, the potential 21 is not determined by and 
a; it has to be calculated independently. Conceptually this is an important implication of 
the underlying quantum geometry. 

We will refer to gab as the dressed, effective geometry to distinguish it from the metric gab 
defined by the effective trajectories, discussed in section II B, that track just the peak of the 
wave function \l/o. Similarly, we will refer to 21 as the dressed, effective, external potential. 
The result on exact equivalence of dynamics holds only for the dressed effective metric 
cjab- However, in practice, from the bounce until the effective trajectory enters the general 
relativistic regime, the wave function \l/o is so sharply peaked on the effective trajectory 
that the difference between gab and gab would be too small to track accurately in numerical 
simulations (since even in sophisticated simulations the numerical errors are much larger). 
Therefore in numerical simulations in section V we will not distinguish between gab and gab- 

Finally, although its coefficients (in the </>, x chart) depend on gab is a smooth tensor 
field. Therefore, it is now straightforward to pass to the conformal time fj it defines. We 
have: 

gabdx^dx^ = ds^ = a{-drf + dx^) , (2.15) 

where 

df, = [e~a']p-^^^d(t>. (2.16) 

This description is especially convenient in the general relativistic regime to make direct 
contact with the cosmology literature. 

To summarize, the dynamics of quantum fields 7^ on a quantum geometry \&o is math- 
ematically the same as that of their dynamics on curved background geometry (2.15). For 
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the scalar mode Qg, the classical 'external potential' 21 has also to be replaced by a quantum 
corrected external potential given by (2.14). This equivalence greatly simplifies our task of 
defining the appropriate Hilbert space and operators for Qg, 7^ because we can import into 
quantum field theory on quantum FLRW space-times of [48] the rich set of techniques that 
have been developed in the QFT on classical FLRW space-times. We will summarize the 
main results in section II C 2. However, we emphasize that the equivalence holds only if the 
perturbations can be regarded as test fields; i.e., their back reaction can be neglected. For 
now we will assume that this is the case. At the end, in section VI, we will carry out a 
self-consistency test by verifying that the assumption is in fact met. 



2. Hilbert space Hi of perturbations and operators thereon 

In this subsection we briefiy summarize the characterization of the physical Hilbert Hi 
of perturbations and the appropriate regularization of relevant composite operators, such as 
energy density. We will work with the adiabatic regularization approach [68, 69], which is 
particularly convenient to perform explicit computations in the cosmological context. For 
simplicity of presentation we will first discuss the tensor modes and then summarize the 
modifications needed to incorporate the scalar mode. 

We can follow the standard procedure used in the cosmology literature since the pertur- 
bations can be regarded as propagating on a FLRW metric Qab-^ The symmetries of this 
dressed, effective background gab allow us to expand the field operator T(x, fj) in Fourier 
modes 

J^si^kek{v) + Al.eUv)) e^'^ , (2.17) 

where, as usual, k = \k\. The operators T{x,fj) satisfy the Heisenberg equation of motion 
if and only if the mode functions efc(r/) are solutions of the wave equation 

el{fj) + 2ie'M + eek{fi) = 0, (2.18) 
a 

were the 'prime' denotes the derivative with respect to fj. These solutions Ckifj) are to 
provide a generalization of the positive frequency basis e"*''*/ \/2k used in Minkowski space- 
time. The canonical commutation relations for the field operator 7'{x,fi) and its conjugate 
momentum imply 

[i,^,it,] = ihi2iT)' 6ik - k') {ckiv), ek'iv))-' and [A^,Ap] = 0, (2.19) 

where _^ 

{eM, ek'iv)) := ^{eM et,{fi) - e'Mel,{fj)) . (2.20) 
Ak 

Therefore, if we require that the basis functions efc(r/) in (2.17) to satisfy the normalization 
condition (6^(77), 6^(77)) = i, then Aj^ and At satisfy the familiar commutation relations of 



^ For a more mathematical treatment involving representations of Weyl algebras, sec [19]. 

^ Since T(x, fj) is purely inhomogcncous in the classical theory, it is natural to require the same property 

for 77). Then, the integral excludes the point fc = 0. But we will require that 6^(77) is continuous in 

k and k 6^,(77) admits a limit as A: —J- 0. 
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creation and annihilation operators. Note that, because 6^(77) and 6^/(77) are solutions of 
(2.18), the scalar product (6^(77), 6^/(77)) is constant in time. Therefore it suffices to impose 
it at some initial instant. 

The first steps in the construction of the Hilbert space "Hi of tensor perturbations can be 
summarized as follows: i) Choose a family of normalized solutions 6^(77) to (2.18); ii) Define 
the associated vacuum state |0) as the state annihilated by all Ag; and, iii) Construct the 
Fock space generated by a repeated action of creation operators At on the vacuum. With 
further important qualifications discussed below, this Fock space will be the required Hi. 

The vacuum state constructed in this way is invariant under the translational and rota- 
tional isometries of Qab- This property will be will be important in our later discussion and 
follows from the symmetries of the 2-point function 

/d'^h 
__ [e^'^'^-'--^ e,{f,{)el{f„)\ . (2.21) 

which suffices to completely characterize this Fock representation. However, it is clear from 
the above construction that the vacuum state is far from being unique: Different choices 
for the bases lead to different definitions of operators A^: and At, and therefore different 
vacua. In the absence of additional physical inputs — such as invariance under the much 
larger isometry groups available in Minkowski or de Sitter space-times — one cannot single 
out a preferred vacuum. More importantly, there is no a priori guarantee that the vacuum 
chosen by a basis would belong to the Fock space determined by a different basis e^. 
There are two sets of potential problems: ultra-violet and infra-red. 

In cosmological backgrounds the ultra-violet issues can be faced by imposing appropriate 
regularity condition on the basis that select the Fock space. In the FLRW backgrounds, 
one can impose the adiabatic condition [68-70], i.e. require that the Ck must approach 
the Minkowski space-time positive frequency modes, e~*'^*/v^2^, at an appropriate rate in 
the limit in which the physical momentum k/d is much larger than the energy 
provided by the space-time curvature. (For a succinct discussion of the technical statements, 
see the companion paper [19].) The rate of approach determines the adiabatic order of the 
basis Ck and controls the ultraviolet behavior of states in the Fock space it selects. Since the 
stress-energy tensor is a dimension 4 operator, for it to be well-defined one needs the basis 
Cfc to be of 4th adiabatic order. We impose this restriction. 

Had the spatial hypersurfaces been compact — with, say, topology rather than M'^ — 
the Fock representations arising from any two 4th order adiabatic bases and e^i. would have 
been unitarily equivalent. In particular, all the 4th order adiabatic vacua would lie in the 
same Hilbert space, which we could take as our "Hi. However even in this case there would 
be no 'preferred' vacuum state: each regular basis would define one and, without additional 
physical input, there is a no preferred basis. However, the total number of particles (0| N |0) 
defined by (the creation and annihilation operators associated with) a regular basis e^ in 
the vacuum state |0) defined by another regular basis Ck would be finite. 

In this paper, the spatial topology is M"^ and there is a further, infrared difficulty. Now, 
the adiabatic regularity only ensures that the spatial number density of the 'under-barred' 
particles in the vacuum |0) is finite. Furthermore, since each vacuum is translationally 
invariant, this number density is constant in space. Therefore, the total number is infinite 
for the trivial reason that the spatial volume is infinite! This infrared infinity is generally 
regarded as physically spurious. However, mathematically it implies that, generically, the 
vacuum |0) defined by a regular basis e^ would have infinite norm in the Fock space defined 
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by another regular basis e^, whence the two Fock representations are unitarily inequivalent. 
However, they can be regarded as physically equivalent in the sense that the expectation 
values of any 'under-barred' observable that refers to a compact region is well-defined in |0) 
(and hence on a dense subspace of the corresponding Fock space). This region could be 
Mlsj the portion of the (infinite) surface of last scattering that is accessible to observations 
(physically, Mls is intertwined with the background quantum geometry \&o)- For example, 
for the number operator N(Mls); corresponding to the number of 'under-barred particles' in 
the spatial region 9^lS) (0| N(D^ls) |0) would be finite for any regular bases Ck and e^.. 

Thus, in the topology under consideration, for definiteness, we will fix a basis which 
is regular up to the 4th adiabatic order and take "Hi to be the Fock space it defines (specific 
examples of convenient bases will be discussed in section IV). All operators of physical 
interest to us will be well-defined on "Hi. Furthermore, observables such as energy density, 
or particle number and energy in a compact region, so constructed, will have well-defined 
expectation values on Fock spaces constructed from any other basis e^ which is also regular 
up to 4th adiabatic order. 

The expansion (2.17) of T(i;, f/) immediately implies that it is a well-defined operator 
valued distribution on "Hi, i.e. that / cPxT{x,f])f{x) is a well-defined (self-adjoint) operator 
on T-Li for every smooth function f{x) of compact support. However, since observables such 
as energy density involve a product of these operator valued distributions, they have to be 
regularized. It is here that the 4th order adiabatic regularity plays a crucial role. It provides 
a natural, mode by mode subtraction scheme that removes the ultra-violet divergences in 
a local and state independent manner, while respecting the covariance of the underlying 
theory. 

Let us begin with the energy density operator for tensor modes. Recall first that in the 
classical theory the stress-energy tensor is given by 

Tab = ^ [VaTVbT- ^^a^r' V,rV,r] . (2.22) 

Therefore, (2.17) implies that the expectation value of the energy density operator p^'^\x, fj) 
in the vacuum state associated with a regular basis Ck is formally given by 

h f d^k 

(0|p(^^|0)formal := (0 1 T.fcr?' 1 0) formal = ^ / 77^H\^ + k^k]^] , (2.23) 

where is the unit normal to the homogeneous slices w.r.t. Qab- The formal expression is 
ultraviolet divergent. The adiabatic regularization scheme provides a specific mode by mode 
substraction, yielding 

(0|p(^)(x,r^)|0) = ^|^[|elP + fc^|e.p-4KC(^)(fc,r^)] , (2.24) 

where the subtraction term C^^\k,fj) is given by 

C^^\k, r^)= + + . (2.25) 



2d^k 8d^k^ 

This subtraction tames the ultraviolet divergences. But there is a further subtlety, again 
because the spatial topology is (rather than T^) and T{x, fj) is effectively a zero rest 
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mass field: the right side of (2.24) has an infrared divergence, directly inherited from the 
term in C^'^\k,fi). Therefore, we need to introduce an infrared cut-off. A natural 
strategy to handle this issue is to absorb modes with wavelengths larger than the radius 
^H(f/) of observable universe into the definition of the homogeneous background itself. Given 
an instant of time fj, this provides us with a minimum physical frequency (or maximum 
physical wavelength) for the modes that are to be treated as arising from perturbations. 
Thus, we are led to the following strategy: Perturbations whose physical wavelength is 
shorter than 9^(r/) will contribute to (O"*^^! p |0°'^^)rcn while those with larger wave lengths 
will be treated as contributing to the background. This procedure is supported by the fact 
that, if we were to use a 3-torus spatial topology, the current observational limits on the 
radius of the torus is close to ^H(t) [71]. Therefore, had we used a torus with this radius, the 
infrared cut-off we use would have been implemented automatically. 

Note that, in terms of co- moving wave numbers k used in our analysis, this infrared 
cut-off is in fact time-independent: It is given by ko = A;^/8.58 where /c^ is the reference 
mode used by the WMAP (see section III A below). As discussed in section VI B of [19], an 
infrared cut-off with fixed co-moving /Cmin leads to a renormalized energy density satisfying 
various criteria that are generally used to select viable renormalization procedures in QFT 
in curved space-times. Thus, the strategy is consistent with the conceptual framework of 
the renormalization theory. 

The final expression of the expectation value of the renormalized energy density operator 
is given by 

(0|p(^)(f,^)|0).en= ^ / [\e'f + e\e,\'-4nC^^\k,v)] . (2.26) 

where C^^{k,fi) is again given by (2.25). Note that the right hand side is independent of 
X because the vacuum is translationally invariant. The same procedure yields the matrix 
elements of p^en on the obvious dense subspace of "Hi obtained by operating on the vacuum 
by finite linear combinations of products of creation operators. These of course depend on 
both X and fj. The matrix elements of the entire stress-energy operator can be constructed 
in the same fashion. We have focused on energy density because, as noted below, it is 
the integrand of the Hamiltonian operator; its expectation values determine whether our 
truncation approximation is valid [19]. 

Remark: As shown in [19], the Hamiltonian operator that governs the evolution of 7^ in 
conformal time fj has the formal expression: 

^(T) _ f d^k ,2k (r),2 Q'feV|2x >.7^ 

^l,fonnal ' Jj^A^ l^k I + ^I^aII ) • (2.27) 

(In light of (2.4), this expression can be anticipated once we know the equivalence between 
quantum fields T propagating on the quantum geometry \E'o and gab-) Now, the standard 
relation between the energy density and the Hamiltonian generating evolution in the 
conformal time implies -^{^ormai — I ^-^^ PkJmai- Therefore, even if we replace p^OTmai 
with prcn with its infrared cut-off /cmim the vacuum expectation value of the Hamiltonian 
operator still has the trivial divergence simply because the integrand is constant and the 
total volume is infinite. However, as discussed above, we can restrict the space integral 
just to the spatial region Mobs(?7) that is observable at time fi and obtain a Hamiltonian 
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operator Hohsiv) tailored to Robs(^)- It provides the desired evolution of operators smeared 
with (space-time) test functions with support anywhere in the observable universe. (For a 
more mathematically complete discussion, see [19].) 

Finally, let us consider the scalar perturbations Q. Conceptually the situation is parallel. 
The mode functions (analogous to the for tensor perturbations) now satisfy the equation 



(2.28) 



The energy density operator (that is needed to test the validity of the truncation approxi- 
mation [19]) is given by: 



2S2 



q'u? + {^ + k') 



k ■ (27r)3 

where the subtraction term C^^\k,fi) is now given by 



C^^\k,fi) 



(2.29) 



C^Q){k,v) = - 



k a'^ + a^^ -a^Qi^ + 2a^a'Qi' + 4a'''h" + a{d"^ -2a'{Qla' + a')) 



2a^k 



+ 



(2.30) 

In summary, we have lifted the adiabatic techniques from quantum field theory in curved 
space-times to quantum fields T and Q propagating on FLRW quantum geometries ^o- 
These techniques led us to quantum states of perturbations that have a good behavior in 
the ultraviolet and to a procedure to systematically regularize products of operator valued 
distributions that are of direct physical interest. Because of M'^ topology, there are some 
subtleties associated with infrared divergences. However, they can be handled by restricting 
attention to spatially compact regions of direct physical interest. 



III. INITIAL CONDITIONS 

Now that we have well-defined Hilbert spaces "Ho and "Hi and physical operators thereon, 
given an initial state on T-Lo ® 'Hi, we can evolve it all the way through the inflationary era 
and calculate power spectra and spectral indices. Our primary goal is to explore if there 
exist viable candidates of initial states that lead to: i) a slow roll inflationary phase that is 
compatible with the 7 year WMAP data, ii) predictions for power spectra and spectral indices 
that are compatible with current observations, and, iii) deviations from the BD vacuum at 
the onset of inflation. If the answer to the first two questions is in the affirmative, we will 
have a viable extension of the inflationary scenario to the Planck regime. If the answer to 
the third question is also in the affirmative, we will have means to test LQC signatures of 
pre-inflationary dynamics for future observations. Since the bounce replaces the big bang 
in LQC we will specify the initial state \Efo ® ^ at the bounce. 

The emphasis in this paper is on investigating the existence of such initial conditions 
rather than on their uniqueness. Nonetheless, at the end, we will provide strong motivation 
for our choices using symmetries, regularity and the novel 'repulsive force' in the Planck 
regime of LQC, created by the underlying quantum geometry. We hope that future investi- 
gations will develop these considerations into a more systematic procedure to arrive at our 
initial conditions from first principles. 
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A. Initial conditions for the background quantum state 

As explained in section II B, in this paper we focus on quantum states of the back- 
ground geometry which are sharply peaked on effective trajectories of LQC Therefore, to 
specify initial conditions for \l/o we have to examine the effective solutions in some detail. 
Our task is two-fold: First isolate the initial data at the bounce for effective solutions, and 
then single out the portion ^ of the permissible initial data set that yields trajectories 
compatible with the 7 year WMAP data. 

Let us begin by noting that in the context of inflationary models, the seven year WMAP 
data [72] is generally parameterized assuming that the state of perturbations at the on- 
set of inflation is the BD vacuum. In this parametrization, a combination of theory and 
observations fixes the infiaton mass to be [6] 

m = 1.21x10"^. (3.1) 

It also implies that the reference co-moving mode used by WMAP® exits the Hubble 
radius during slow roll at a time at which the Hubble parameter H, the slow roll 

parameter e = —H/H"^, the infiaton 0, and its time derivatives (p have the following values 
[6]: 

H{t{h)) = 7.83 X 10"^; and e{t{h)) = 8 x 10^^ 

(f)it{h)) = ±3.15 and ^{t{h)) = 1.98 x 10"^. (3.2) 

where the 'dot' refers to the cosmic time. Because of observational error bars, these values 
are uncertain in a ~ 2% window. We will use these values in our analysis and ask at the end 
whether the quantum state of perturbations at the onset of the slow roll is indistinguishable 
from the BD vacuum. In most of the parameter space that dictates the initial conditions at 
the bounce, the answer will be in the affirmative. For a narrow window in which it is in the 
negative, one would have to re-calculate values of these parameters (see Section VA). 

Our task then is to find the class of initial conditions for the effective equations for which 
dynamical trajectories enter the tiny region of the phase space defined by (3.2) and the 
associated error bars, some time in their future evolution. As in general relativity, the space 
of initial data in the effective theory is 4 dimensional: specification of z/, b; 0, at any time 
determines a unique solution to effective equations. (Recall that z/ ~ a'^ and H ~ sin2Ab.) 
At the bounce, we have Ab = 7r/2 and furthermore, as in general relativity, we have two 
symmetries on the space of solutions: i) a rescaling symmetry , z/(t) — > az/(t), (f){t) — )■ (f){t), 
and, ii) and the time reflection symmetry vit) i>{—t), (j){t) — (f){—t). The first allows us 
to restrict ourselves to solutions with z/|b = for some fixed constant u", while the second 
lets us focus just on solutions with 0|b > 0; solutions with arbitrary i/|b and also those with 
0|b < can be obtained using these symmetries. Finally, the LQC Friedmann equation 
(2.10) implies that the value 0b of at the bounce point determines |0|b there. Thus, the 
free data at the bounce is just the value 0b of the background infiaton. Next, since the total 



fc* is given by /c^/co = 2 x lO^'^Mpc"^, or, fc^ = 8.58 fco, where Uo refers to the scale factor today and, 
as before, ko refers to the wave number that has just re-entered the Hubble radius today. It is only the 
combination fc*/ao that has direct physical meaning; 2T:ao/ki, is the physical wave length of this reference 
mode today. 
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energy density at the bounce is fixed, p = pmax on all solutions, it follows that this free data 
(pB is restricted to lie in a finite interval: 

|0b| e (O,0max), wherC 0max = ^ ^ 7.47 X 10^ (3.3) 

m 

This interval provides us with the relevant parameter space in this paper. 

For definiteness, we will assume 0b > since the sign does not make a qualitative 
difference to the analysis. In this case, it was shown in [6] that a sufficient condition for the 
trajectory to enter the tiny phase space region compatible with the 'WMAP slow roll' is that 
(pB > 0.93. This excludes only a tiny part of the full range of allowed values. In this sense, 
in the effective theory, 'almost all' initial conditions at the bounce point are compatible 
with the 7 year WMAP data. We would like to emphasize that this behavior is very non- 
trivial. At first sight it may seem that the statement just says that, as in general relativity, 
infiationary trajectories are attractors in effective LQC. However, the implication is much 
stronger because compatibility with the seven year WMAP data is a sharp quantitative 
requirement: We need the trajectory to achieve (p = 3.15 and e = 8 x 10"'^ at the time when 
H = 7.83 X 10~^ within the WMAP error bars. Second, as we saw in section II B 2, pre- 
inflationary dynamics of effective LQC have several features that are very distinct from what 
occurs in general relativity. Therefore, a priori one cannot assume that even the attractor 
property of general relativity must carry over. Finally the pre-inflationary dynamics covers 
some 11 orders of magnitude in matter density and curvature. Therefore it is rather striking 
that 'almost all' effective trajectories starting at the bounce would enter the tiny region of 
phase space compatible with the onset of the desired slow roll. 

From now on we will restrict ourselves to this part, 0b ^ 0.93, of the parameter space. 
The restriction involved is tiny. However, this does not imply that the initial conditions for 
the quantum state \l/o at the bounce are generic. The permissible quantum states are very 
special because we have asked that \l/o be sharply peaked at a point on the constraint surface 
of the phase space of the effective theory at the bounce time. What is generic is only the 
point at which they are peaked. 

We will conclude by clarifying a subtlety about what it means to have the initial state 
^o(i^, 0b) at the bounce to be peaked at a point (i^°,b = vr/2A; (Pb^V"^)) the constraint 
surface. Note first that for any given 0g and z/°, the Hamiltonian constraint (2.10) in 
the effective theory determines The subtlety is that, as noted in section II B, the 

Hamiltonian constraint is de-parameterized on the physical sector: Since serves as 'internal 
time', it is a parameter rather than an operator on the physical Hilbert space T-io- However, 
because of the quantum Hamiltonian constraint, P{(i>)^o = —Ho^oi Piip) is a well-defined 
operator on "Hq. Therefore, the initial state \l/o(z^, 0b) is to be chosen so that the expectation 
value of v is that of Ho is the determined by the chosen 0g, and fiuctuations in both 
quantities are very small. For further details on construction of such states, see [22] . 

B. Initial conditions for the quantum state if) of perturbations 

For linear test fields in de Sitter space-time, one can single out a unique state by demand- 
ing that it be regular and invariant under the full de Sitter group. This is the BD vacuum. 
Since in the slow roll infiationary phase the Hubble parameter is approximately constant, 
space-time geometry is well approximated by the de Sitter metric. This motivates the usual 
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choice of the BD vacuum as the quantum state of perturbations at the onset of inflation. 
However, in the pre-inflationary phase, especially near the bounce, space-time geometry is 
very far from the de Sitter geometry. As we noted in section II B 2, the Hubble parameter 
vanishes at the bounce. It then increases very quickly to attain a Planck scale value (in less 
than a Planck second for kinematic energy dominated bounces!) [6]. Therefore, it is not 
possible to even formulate the notion of an approximate BD vacuum at or near the bounce. 

However, we can still carry over the central idea that led to the choice of the BD vacuum in 
the standard inflationary scenario and ask that the initial state Q%) of perturbations 

be regular and maximally symmetric. Since the state o and hence the dressed effective 
metric (jab is invariant under the FLRW isometries, the symmetry requirement is precisely 
that V'('7i; Qfc) be invariant under spatial translations and rotations. We are thus led to 
ask that 'ip{Tj:, Qg) be a 4th order adiabatic vacuum at the bounce time fj = fj-Q. In the 
Schrodinger picture, the state will evolve but retain its 4th order regularity and invariance 
under rotations and translations (see Eq (2.21)). In the Heisenberg picture, of course, it is 
manifest that the conditions are not tied to any specific time. 

Although we have imposed maximal possible symmetry requirements that are compatible 
with pre-inflationary physics, as our discussion of section II C shows, the space of permissible 
states is still infinite dimensional. A second requirement on the choice of the initial if) is 
motivated by our truncation approximation. Clearly, for the approximation to be meaning- 
ful, the initial perturbed state ip should be such that its contribution (t/'I p(x, ?7b) |'^)ren to 
the energy density at the bounce time f/B should be negligible compared to Pmax; the energy 
in the background quantum geometry In contrast to the regularity and symmetry, this 
last requirement is tied to a specific time, fj = fj-Q. 

A key question is whether states satisfying both these conditions exist. Detailed calcu- 
lations using 4th order adiabatic regularity show that they do. A particularly convenient 
choice for numerical simulations can be arrived at as follows. For definiteness, let us con- 
sider tensor modes. Then, in the large k limit, any 4th adiabatic order basis ek{fi) has to 

approach the explicitly known approximate (WKB) solutions fj^^\'fj) to (2.18) faster than 



(S/fc)9/2. The /i'^(f/) are given by [19] 



fl'\v) = , e-/'^^i''(-)'^- (3.4) 



where wi^\fj) = Eti'^i' with 



1 a' 



Wo = k; Wi = 0; W2 = - — -^ 

2k a 



2a" a'"^ — 2a"^a — 2aa'a + a^a 
W:, = 0; W,^ (3.5) 

The leading order term of (3.4) corresponds to the positive frequency solution in Minkowski 
space and the rest of the terms are higher order contributions that vanish at different rates 
when (a/k) — )■ 0. Therefore, a natural strategy is to construct an 'obvious' 4th order regular 
basis e^^^(f/), tailored to the bounce time fj^, as follows. Although the definition of 4th order 
adiabaticity is only an asymptotic one, we can ask that ef'^{fj) and f^f'i'f]) share the same 
initial data at the bounce time fj = fj-Q: 

ef^(r?B) = flt\m) ; and a^ef^(r/B) = d^fP{f,^) . (3.6) 
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(This relation will not hold at future times because e^^^ satisfies (2.18) exactly while fj, 
satisfies it only approximately.) The vacuum state |0°^^) determined by e^^^ will be referred 
to as the 'obvious 4th adiabatic order vacuum'. 

There is one subtlety: For this construction to work, W^'^\fjB) in (3.4) must be non- 
negative since it appears under a square-root in the expression of /^^^(r/e)- This is the case 
when the infra-red cut-off is not too low. In this case, detailed numerical calculations have 
shown that (0°^^| p(x, r^e) |0°^^)ren is less than a few percent of the background energy density 
Pmax- If the infrared cut-off is lower, on the other hand, one has to modify the procedure 
to obtain the initial data for e^^^: It is to be given by a suitable smooth extension of /^'^^ 

and its time derivative from the high k values where M^/^^^(?7b) is non-negative to the lower 
k values where it becomes negative. Consequently, jO"*^^) is not as 'canonical'. However, 
the required extension can again be carried out keeping the renormalized energy density in 
perturbations only a few percent of that in the background. 

The non-trivial question is whether the energy density in the perturbation continues 
to remain negligible for fj > fj-Q. The issue is especially non-trivial in the Planck regime. 
However, using numerical evolution we will show in section V that the answer is in the 
affirmative for the state |0°^^) at least when 0b > 1.23. Furthermore, we will argue 
analytically that the same is true for states in an open (infinite dimensional) neighborhood 
of |0°*^^). Thus, states satisfying all our initial conditions at the bounce do exist. 

Remark: Physically, the IR cut-off is dictated by the physical radius 9^ls of the observable 
universe at the surface of last scattering (see section II C 2). More precisely it is the co- 
moving radius 9^ir such that £Hls = (^^(^ls)) 5^ir- Since £Hls is observationally fixed, the 
value of 5^iR is sensitive to the number of e-foldings between the bounce and the surface of 
last scattering. Now, since the background states ^ o under consideration all meet conditions 
that are compatible with the 7 year WMAP data at the onset of slow roll, the number of 
e-foldings between this onset and the surface of last scattering is the same for them. What 
differs from one o to another is the number of e-foldings between the bounce and the onset 
of slow roll; they grow rapidly as 0b increases. Therefore the subtlety in the choice of mode 
functions mentioned above arises only for low values, 0b < 1.18, for which the number of 
e-foldings between the bounce and the onset of inflation is low. 



C. Physical consideration 

In this paper, our primary goal is to explore the extent to which the inflationary scenario 
can be extended to the Planck regime. We saw that, in effective LQC, a very large portion 
^ of the allowed initial data at the bounce leads to dynamical trajectories that encounter a 
slow roll phase compatible with the 7 year WMAP data. This led us to choose an initial state 
'^o of the background geometry to be sharply peaked at a point of The initial conditions 
on the states ip of perturbations were motivated by three considerations: i) Regularity; 
ii) Symmetry; and iii) Compatibility with our truncation approximation. The flrst two 
conditions led us to restrict V to be a 4th adiabatic order 'vacuum' state and the third 
led us require that the energy density in if) be negligible compared to that in background 
quantum geometry at the bounce. This set of conditions is not overly restrictive in the 
sense that they allow an infinite set of initial states. At the same time, because '^o is required 
to be sharply peaked at a point in and because ip has to satisfy the three conditions at 
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the bounce, the allowed initial states '^o ^ constitute a very tiny subset of all states in 
"Ho ® "Hi- In this sense they are very special. We will now provide heuristic considerations 
to help clarify the physical motivation that underlies our choices. This discussion could also 
serve to further narrow down the initial conditions and opens some avenues for further work. 

A common strategy, articulated by Penrose [73] in particular, is to use known physics and 
our knowledge of the current universe to draw conclusions about the special nature of initial 
conditions. In LQC this general idea has been used to narrow down the choice of the initial 
quantum state of the background geometry [7]. As noted in section III A, one begins 
with the observation that, away from the deep Planck regime, the universe is extremely 
well- described by a FLRW solution of Einstein's equation. One therefore takes states in Tio 
which are sharply peaked on a FLRW solution and evolves them back in time. One finds 
that the states continue to be sharply peaked with small uncertainties in both 'conjugate' 
variables Ho and 0, all the way to the bounce. It is this property that motivates our initial 
conditions on \E'o at the bounce. 

The choice of tp can be further narrowed down by similar considerations. In any of our 
vacua |0) for perturbations, the expectation values of 7^, and their conjugate momenta 
vanish, encoding the idea that departures from homogeneity and isotropy are small. But we 
could also ask that the states be 'as peaked on the zero values as is allowed by the uncertainty 
principle'. There is a precise sense in which the state 10°*^^) we discussed in some detail in 
section IIIB satisfies this condition. Heuristically, this condition would say that the initial 
state ® is as homogeneous and isotropic as is permissible in the quantum physics of 
the truncated theory. Therefore this condition will be referred to as quantum homogeneity 
and isotropy at the bounce. A precise formulation and implications of this requirement will 
be discussed in detail in [74]. 

As we explained in the beginning of this section, our goal is only to show the existence 
of initial conditions that lead to a self-consistent completion of the standard inflationary 
paradigm and we will see in the next two sections that quantum homogeneity and isotropy 
at the bounce meets this goal. In future and more comprehensive investigations one could 
inquire if one can arrive at this condition from deeper considerations. Specifically, since 
the background quantum geometry \E'o has a pre big bang branch, one can ask if the initial 
conditions at the bounce naturally arise from the prior evolution. We will conclude by 
providing some heuristics that suggest a promising direction to address this issue. 

First, if 03 > 1.15 (the value used in most numerical plots in this paper), the universe that 
is observable at the surface of last scattering expanded from a ball of radius re < 10£pi at the 
bounce. Therefore, to account for the CMB observations, one needs quantum homogeneity 
and isotropy at the bounce only at this small scale. At first, this smallness of scale may 
appear to make the requirement quite weak. Indeed, such an argument is sometimes made 
in the standard inflationary scenario where one assumes the validity of general relativity all 
the way to the Planck scale and argues that a similar small region expands out to fill the 
universe that is observable at the surface of last scattering. However, it has been pointed 
out by many authors, on the initial ball the required homogeneity and isotropy has to be 
extraordinarily high and a natural mechanism to achieve this does not appear in the standard 
scenario. We will now give qualitative arguments to the effect that the underlying quantum 
geometry of LQG may provide the missing element. 

Prior to the bounce, the universe would be collapsing thereby producing strong inhomo- 
geneities that lead to complicated growth of curvature. However, the LQC models that have 
been analyzed in detail suggest the following overall picture. Quantum geometry effects 
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lead to a new repulsive force. This force is completely negligible until a curvature scalar 
approaches the Planck regime. But then it grows very quickly, overwhelms the classical 
gravitational attraction that would have made the curvature scalar singular, and dilutes its 
value. Thus, in addition to the 'global' density bounce we have focused on, there would be 
many 'local' bounces associated with local growth of curvature. The quantum geometry cor- 
rections start to be significant when the curvature is ~ 10~^£pf, (i.e. about a thousandth of 
the maximum value) which corresponds to a radius of curvature of about 10 ipi. Therefore, 
one may expect the 'dilution' mechanism to be effective in ironing out the 'wrinkles' in the 
curvature also on that scale. Thus, the missing element in achieving the required quantum 
homogeneity and isotropy on the 10£pi scale could well be provided by the repulsive forces 
that originate in quantum geometry. (See [64, 65] for a semi-classical mechanism leading 
to desired initial conditions in a different context.) This mechanism would not preclude 
inhomogeneities and anisotropies on larger scales at the bounce. But these would corre- 
spond to modes whose wavelength at the surface of last scattering is much greater than the 
radius of the observable universe and would not be observable. If these considerations are 
borne out, the pre-bounce history of the background would not have a direct relevance to 
the post-bounce evolution of the portion of the universe that is observable. 

Of course a much more detailed investigation of this 'dilution effect' is needed to deter- 
mine if these qualitative ideas are viable and can be developed into a detailed quantitative 
argument. 

Remark: If we just focus on the tensor perturbations, then the requirement of quantum 
homogeneity and isotropy at the bounce implies that the state ^'o(^) ® satisfies a 

quantum version of Penrose's Weyl curvature hypothesis [73] in the following sense. In its 
original formulation, space-time geometry was assumed to be classical, and the condition 
was that the initial singularity was very special in that the Weyl curvature vanishes there. 
In LQC of course there is no singularity and hence it is natural to impose a condition at 
the bounce surface at which the matter density achieves its maximum value. Furthermore, 
since the electric and the magnetic parts of the Weyl tensor, Eab, Bab, do not commute 
in the truncated theory, and their commutator is a c-number, there is no state on which 
they can both vanish. Thus, the Heisenberg uncertainty principle tells us that it is not 
meaningful to ask that the Weyl tensor operator be zero even at an instant of time. This 
is completely analogous to what happens in the Maxwell theory in Minkowski space, where 
there is no state at all in the Fock space that is annihilated by the Maxwell field operator 
Fab, or equivalently, the electric and magnetic field operators E"', B°^. The expectation 
values of both these operators can vanish but this can be achieved on a very large class of 
states. But we can impose, in addition, the next natural requirement: i) the product of 
uncertainties in the electric and magnetic fields be saturated; and, ii) the uncertainties be 
divided equally between the two fields. Some care is required in formulating these conditions 
in a precise manner. When this is done, one finds that there is exactly one state that 
satisfies the condition: The standard Maxwell vacuum! In the truncated theory now under 
considerations, it is natural to elevate Penrose's Weyl curvature hypothesis to the quantum 
theory by asking for states ip in which the expectation values of Eab, Bab vanish, the product 
of uncertainties is minimized and the uncertainties are equally divided between the electric 
and magnetic parts. This quantum version of the Weyl curvature hypothesis is satisfied in 
a weaker sense (i.e., only in the adiabatic limit) by states '^o®'4' which meet our quantum 
homogeneity and isotropy requirement. (For details, see [74].) 
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IV. PRE-INFLATIONARY DYNAMICS OF THE QUANTUM CORRECTED 

BACKGROUND 



This section is divided into two parts. In the first we discuss the interplay between the 
dynamics of the background curvature and that of modes of perturbations in physical terms. 
In the second we describe results of numerical simulations of the evolution of the effective 
background geometry. 

To discuss dynamics in detail, one has to choose a time parameter defined by the effective 
background metric. In the cosmology literature one generally uses conformal time rj which 
moves the big bang to 77 = —00. In cosmic time t, it occurs at t = 0. Now, we would like to 
start our evolution at the LQC big bounce which replaces the big bang of general relativity. 
To compare the results with those of general relativity, it is technically more convenient to 
use the cosmic time t defined by the gab and set the bounce time to be = 0. Therefore, 
while our analytical considerations use conformal time fj, our numerical evolutions are all 
carried out in cosmic time t. 



A. Why pre-inflationary dynamics matter 
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FIG. 1: Time evolution of the curvature radius (red solid line with dashes) and of a wave length of 
interest to observations (green dashed line). Left Panel: General relativity. The modes of interest 
have wave lengths less than the curvature radius all the way from the big bang (t = 0) until after 
the onset of slow roll, shown in the inset (t ~ lO^tpi). Right panel: LQC. The bounce occurs at 
time t = and we have set a|j=o 
mode whose wave length Aphy \{ 
modes with smaller wavelengths remain within the curvature radius until the onset of inflation. On 
the other hand, modes with physical wave length greater than the curvature radius at the bounce 
(for example green dashed line) enters the curvature radius soon after the bounce and remain 
within the curvature radius until after the onset of inflation, shown in the inset (t ~ lO^tpi). These 
modes can be excited due to curvature while their wave length is greater than the curvature radius 
and will not be in the BD vacuum at the onset of inflation. (There are two points (at t ~ O.Stpi 
and t w 5 X lO^ipi) in the LQC evolution at which the u; = 1/3 in the effective equation of state, 
whence the scalar curvature vanishes and the radius of curvature goes to infinity.) 
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There is a common lore that inflation would simply dilute away all the effects of pre- 
inflationary dynamics. This belief stems from the following physical considerations. For 
definiteness let us consider tensor perturbations. To bring out the physics that governs 
evolution, it is convenient to simplify the dynamical equation by noting the following 
rnathematical fact. Under a general conformal transformation cjab = VL~^gab we obtain 
(□ - (1/6)^)0 = fi-3(n - (l/6)i?)</) where = fi-^- Now, if we let = a, then gab is the 
fiat metric defined by the coordinates r/, x and we have Dcfc = [□ + (1/6) a?E] (a Ck) = 0. 
Therefore, the fact that 6^ satisfies the wave equation Defc = w.r.t. Qab (see (2.18)) implies 
that the rescaled function Xkiv) = o,iji)^kiji)i satisfies a_wave equation w.r.t. the flat metric 
Qab in presence of a time dependent potential (1/6) a? R: 



effect of curvature on their evolution can be ignored. Put differently, if the physical wave 
length Aphy = 271 a /k of the mode is much smaller than the radius of curvature (6/i?)^/^, 
the mode propagates as if it is in fiat space-time (defined by fjab)] its dynamics is trivial. In 
the standard inflationary scenario, modes that can be observed in the CMB have physical 
wave length smaller than the curvature radius at the onset of inflation.^ If one were to use 
general relativity to evolve back in time, then this inequality can continue to hold all the 
way to the big bang (see the left panel in Fig. 1). Therefore, the modes would remain in the 
same quantum state throughout the pre-infiationary era, whence one would conclude that 
the pre-infiationary dynamics does not have any effect on modes that are observable in the 
CMB. 

What is the situation in LQC? Because of the exact equivalence between the evolution 
of perturbations T on the quantum geometry ^ „ and on the dressed-effective metric gab 
given by Eq. (2.12), these general considerations continue to hold all the way back to the 
bounce. However, the pre-infiationary dynamics of the scale factor a{fi) and the curvature 
R{fi) is now qualitatively different. The key question then is: Do the physical wavelengths 
Aphy of the relevant modes ever become comparable to (6/i?)^/^ during the pre-infiationary 
evolution? If they do, their dynamics would be non-trivial. 

A detailed analysis shows that this is the case for modes whose physical wave length 
Aphy exceeds the curvature radius at the bounce time (see the right panel of Fig. 1). As 
Parker [68, 70] showed already in the sixties, modes that experience curvature are excited. 
Indeed, during slow roll, this is the phenomenon that leads to the inhomogeneities that 
seed the large scale structure. Qualitatively the situation is similar in the pre-inflationary 
epoch as well: modes with Aphy \t=o> iG/RY^'^\t=o are excited in the Planck regime that 
immediately follows the bounce. As a consequence, at the onset of inflation, the quantum 
state of perturbations is populated by excitations of these modes over the BD vacuum. 
Thus, pre-inflationary dynamics can change the initial conditions for perturbations at the 




(4.1) 



where R{ri) is the scalar curvature. Eq. (4.1) brings out the fact that for modes with 
physical wave numbers k/a much larger than the curvature energy- scale kR = {R/Qf/^, the 



For (quasi-) dc Sitter space-times this is equivalent to requiring A <C where dKu = 1/H is the Hubble 
radius. This is the condition generally discussed in the inflationary literature. However, away from slow 
roll, is not so simply related to the radius defined by the curvature scale. 
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onset of the slow roll. Can this lead to physically distinct consequences? The answer is in 
the affirmative. It has been shown that the predictions for the CMB and for the distribution 
of galaxies are sensitive to the quantum state of perturbations at the beginning of inflation 
[10-14]. Furthermore, concrete observational tests to probe these consequences have recently 
been proposed [15-17]. 

To summarize, pre-inflationary dynamics can lead to predictions that are distinct from 
that of the standard inflation. In principle, there may be excitations over the BD vacuum 
at the onset of inflation which are so large as to be in conflict with the observed power 
spectrum and spectral index of the scalar perturbations. If this occurs, the corresponding 
initial conditions at the bounce would be ruled out observationally. If the predictions are 
compatible with current observations, the departure from the BD vacuum can still lead to 
effects that could be observed in the future. This possibility is of special interest because 
it relates the Planck scale dynamics to observations. A priori we do not know if any or 
all of these possibilities would be realized for the permissible values of the free parameter 
of the theory, (ps- These questions can be answered only by a detailed analysis of the 
pre-inflationary dynamics. 

B. Quantum corrected background geometry 

Recall from section II C 1 that, because \l/o is so sharply peaked, within numerical errors, 
one can replace the dressed effective metric cjab with the cjab that traces the evolution of the 
peak of \l/o and satisfies effective equations. This is a set of three equations. The ffist is the 
Hamiltonian constraint (2.11) that must be satisfied at any instant of time and the other 
two provide evolution: 

V{<P)Y + Vmpra,^-V{<P))] (4.2) 
0, 

where as usual the 'dot' denotes derivative with respect to comic time i. As discussed in 
section III A, at the bounce time the only free parameter is the value 0b of the background 
inflaton which is constrained to lie in the interval 0b G (0.93, 7.41 x 10^). For each choice 
of 0B in this interval, we obtain a solution to the effective equations which enters the slow 
roll compatible with the 7 year WMAP data some time in its future evolution. 

The dynamical equations (4.2) constitute a set of coupled ordinary differential equations. 
They were solved using Mathematica. As an accuracy check, the satisfaction of the Hamil- 
tonian constraint (2.11) was monitored throughout the evolution. For reasons mentioned in 
section III A and will become clearer in this subsection, in this discussion we will pay special 
attention to the case when the bounce is kinetic energy dominated. 

The bounce occurs at t = and we use the convention a|f=o — ^ (rather than Stoday = !)• 
Immediately after the bounce there is a phase of super-inflation — i.e. faster than exponential 
expansion — because if > at the bounce and H > immediately after the bounce in 
all solutions. At the end of this super-inflation phase, the Hubble parameter achieves its 
maximum value i^max ~ 0.93. For bounces in which the kinetic energy dominates, the super- 
inflation phase is dominated by the Planck scale dynamics because the Hubble parameter 
grows from zero to its maximum value in a fraction of a Planck second. This extremely 
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FIG. 2: Evolution of the scale factor a (left) and the energy density p (right) as a function of 
cosmic time for (p^ = 1-15. Slow roll occurs around i ~ 10^ tpi and during this phase the scale 
factor increases exponentially and the density remains nearly constant. The insets in both plots 
highlight the behavior near t = 0. The lower curve in the left inset and the upper in the right inset 
show the behavior in general relativity in which the scale factor goes to zero and the energy density 
diverges. By contrast in LQC (upper curve in the left inset and the lower curve in the right inset) 
both quantities remain finite at the bounce. 



short lived phase is followed by a much longer phase during which the LQC corrections 
to dynamics become progressively weaker as the infiaton loses kinetic energy because the 
friction term (due to the Hubble parameter) is large; it lasts till t ~ lO^tpi. At the end 
of this phase, the kinetic energy in the infiaton is equal to the potential energy and the 
total energy density is about IQ-i^ ppi (and hence the curvature has also fallen by ten orders 
of magnitude). The infiaton is still going up the potential at the end of this phase. It 
continues to lose kinetic energy and at t ~ 10^ tpi, it stops climbing, turns around and 
starts descending the potential. Soon there after, when t ~ 10^ — lO^tpi, it enters the slow 
roll inflation compatible with the 7 year WMAP data. The slow roll phase lasts for about 
the same time interval as it takes to reach that phase starting from the bounce through 
pre-inflationary dynamics. Fig. 2 illustrates the dynamical evolution of the scale factor and 
energy density using (ps = 1.15. (See also [6].) 

To make contact with the WMAP phenomenology, we need to locate the reference mode 
ki,. We fix a value of (ps and obtain the corresponding solution a{t), (f){t). We then monitor 
values of the Hubble parameter H{t) and locate the time t{ki,) at which it assumes the 
reference value H{i{K)) = 7.83 x IQ-^ determined by the WMAP data (see (3.2)). Then 
is the co-moving wave number of the mode whose physical wave number equals the Hubble 
parameter at this time: k^, / d{t{k-^)) = 7.83 x IQ-^. This procedure determines both /c^ and 
the time t{k^) at which the reference mode /c^ exits the Hubble radius for any given value 
of (j)B- As Fig. 3 shows, k^, increases very rapidly with 0b • We also calculate e(t(/i;*)) and 
verify that it is in a small window around 8 x 10"^ in order to ensure that the slow roll phase 
is compatible with the 7 year WMAP data within its error bars. As a check on numerics, we 
also verify that (j)(t{k^,)) and (j){t{k^,)) are in agreement with (3.2). Finally, we calculate the 
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FIG. 3: The value of k^, as a function of 0b- As the inset shows, k^, always increases with (p^ but 
the increase is especially pronounced for (p^ > l-l- 
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1.36 X 10~3 


1.97 X 10^ 


20.2 


1.5 


2.7 X lO'^ 


2.3 X lO"'^ 


2.9 X 10^ 


28.9 



TABLE I: This table shows the value of the reference co-moving momentum fc^ used in the WMAP 
data, the corresponding physical wavelength A*(tB) at the bounce, the time t{k-^) at which the 
mode k^, exits the Hubble radius during inflation, and ln[d{t{k^) / d{tB)], the number of e-folds of 
expansion between the bounce and t{k^). We focus on the range for (j)-Q that is relevant to explore 
whether pre-inflationary dynamics can lead to deviations from the BD vacuum at the onset of the 
slow roll. 

physical wave length \^, of the mode ki, at the bounce point and the number of e-foldings 
between the bounce, t = and t = t{k^,). Results for an interesting range of values of 0b 
are collected in Table 1. 

Finally, we can use Table 1 to analyze the effect of the pre-inflationary dynamics on the 
initial state of perturbations at the onset of inflation. The argument is semi-heuristic but 
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its main conclusions provide reliable guidelines because it is based on rather simple physical 
considerations. At the bounce, the value of the background scalar curvature is universal, 
Rb = 62. This provides a new energy scale /clqc and the associated wave-length Alqc : 



/62\ ~ 2n 
kLQC = hr r ~ 3.21, and Alqc = ~ 0.977 , (4.3) 

V / fcLQC 



(see Eq. (4.1) for the origin of the factor 6 and e.g. [75] for a discussion of the observational 
consequences of states with preferred scales). Since a|B = 1, /j^lqc and Alqc are physical 
energy and length scales at the bounce time. Thus, modes whose physical wave-length is 
larger than Alqc at the bounce will experience curvature sometime during their evolution 
and, as Fig. 1 shows, this will occur in a short time interval soon after the bounce. Will 
any of these modes be in the observable range for CMB? Recall that, at the surface of last 
scattering, the physical wave length of the reference mode fc* is about (1/8.58) x ^Hls 
where SHls is the radius of the observable universe at that time. Therefore, we will have 
modes which are both observable in the CMB, and which will experience curvature during 
their pre-inflationary evolution provided the value of 0b is such that 

A, > ^ 0.11. (4.4) 

- 8.58 ^ ^ 

Table 1 shows that this is possible only if 0b < 1.2. If 0b ^ 1-2, none of the modes that are 
in the observational range will encounter significant curvature during their pre-infiationary 
evolution. Therefore, the initial 4th order adiabatic vacuum — say, |0°*^^) discussed in section 
lllB, for definiteness — we start out with, will not be excited in these modes. Therefore 
the state |0°'^^) we constructed at the bounce will be indistinguishable from the obvious 4th 
adiabatic order vacuum constructed at the onset of inflation. But since the physical wave 
lengths of the observable modes are all within the Hubble radius at the onset of slow roll, 
it turns out that the BD vacuum is observationally indistinguishable from the obvious 4th 
adiabatic order vacuum defined at the onset of inflation. Thus, the net result is that 10°*^^) 
we started out at the bounce will have negligible excitations over the BD vacuum |BD) at 
the onset of inflation in the observable modes. If, on the other hand, 0b < 1.2, the state at 
the onset of inflation would carry excitation over |BD) in modes that are observable in the 
CMB and this can serve as a source of new effects [12, 15-17]. This is why we have focused 
on kinetic energy dominated bounces. 

Of course this conclusion can only be taken as a general guideline. Detailed calculations 
are necessary to verify that they did not miss a subtle but important point and to sharpen the 
conclusions through precise quantitative results. In particular, there is no a priori guarantee 
that the power spectrum will be compatible with the WMAP observations for any value of 
0B. Nor is it clear that the underlying truncation scheme is self-consistent, i.e., that the 
back reaction of the quantum perturbations can indeed be neglected all the way from the 
bounce till the onset of slow roll. Section V addresses the first issue, and section VI, the 
second. 



V. PRE-INFLATIONARY DYNAMICS OF PERTURBATIONS 

In this section we will summarize the results of numerical evolution of quantum fields 
representing scalar and tensor perturbations all the way from the bounce to the end of slow 
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roll inflation. We will then calculate the scalar and tensor power spectra and compare them 
with those of the standard inflationary scenario. This discussion will provide a conceptual 
completion of the scenario to include the Planck regime. In addition, we will flnd that the 
pre-inflationary dynamics can provide corrections to the standard inflationary scenario with 
potentially observable consequences. 

This section is organized as follows. In VA we flrst recall the WMAP phenomenology 
and then summarize the modiflcations that arise if the state is not the BD vacuum. In V B 
we sketch the plan of the calculations. Illustrative choices of initial conditions are discussed 
in VC and the results of numerical evolution are presented in VD. Section VE discusses 
a sample of checks that were performed on the numerics and VF summarizes the overall 
picture. In broad terms, detailed numerical calculations bear out the physical scenario 
developed in section IV B and provide detailed quantitative information on the scalar and 
tensor power spectra in LQC and their relation to those in the standard inflation. 



A. WMAP: interplay between theory and observations 

As noted in section III A, the parametrization of the WMAP data is normally carried out 
assuming that the state of perturbations at the onset of inflation is the BD vacuum. We will 
start with this parametrization. For values of 0b that lead to the BD vacuum at the onset, 
as in the standard theory, we will have a self-consistent analysis. When there are signiflcant 
departures from the BD vacuum, one has to revisit the issue and flnd a parametrization that 
is consistent with the new quantum state at the onset. 

Let us begin with the scalar perturbations. To compare with observations, one gener- 
ally calculates the power spectrum V-jiiM) of co-moving curvature perturbations TZ^ (see 
Appendix A). This is because they are conserved to an excellent approximation once the 
mode exits the Hubble radius during the slow roll: 7?.'^ — )■ for k/a < H. This insensitivity 
to the details of the post-inflation dynamics of the background geometry greatly facilitates 
the task of relating the spectrum of perturbations at the end of inflation with the CMB 
temperature fluctuations. It is therefore tempting to carry out the entire analysis starting 
from the bounce using TZj^. Unfortunately, these flrst order gauge invariant perturbations 
fail to be well-deflned at the 'turn around' during the inflaton evolution where = 0. On 
the other hand, the Mukhanov- Sasaki variables Qj: introduced in section II are well deflned 
throughout the evolution from the bounce to the end of inflation. Therefore, we numerically 
evolve Qj: starting from the bounce, calculate the power spectrum VQ{k), and convert it to 
the power spectrum V-jzik) at the end of inflation via 



z 



a 



2 



6 P{^) _ 7 



VQ{k) = - Vn{k), where z{f,) = — £11^ ^ £1^ . (5.1) 



K 7r(a) a?P b 



Given a 4th adiabatic order vacuum of the scalar perturbation, the power spectrum 
VQ{k) for the Mukhanov- Sasaki variables is deflned in terms of the associated two point 
function in the momentum space: 



{n\ Q^.Qp \n) := {2ii)H{k + k') ^V^{k) , (5.2) 



where all quantities are evaluated a few e-folds after the time 77^, deflned by k/a{r]k) = H{r]iS), 
at which the mode k 'exits the Hubble radius'. As discussed above, at the onset of the slow 
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roll we will first set \Q) = |BD), the BD vacuum, and introduce the modifications that are 
needed for a general state \Q) at the end. As is well-known, a basis qf^ifj) that defines |BD) 
is given by 



-krj) 



(5.3) 



where Hj^\—kri) is a Hankel function (of the first kind). The index n is determined by the 
slow roll parameters, fi = 3/2 + 2e + 6, with 



H 



2HH 



1 ,V_, 



V^2 V ' ^ V 



(5.4) 



Here the 'dot' refers to the derivative with respect to the cosmic time t and we have 
ignored quantities quadratic in slow roll parameters to approximate them by functions 
of the potential V{(j)). The mode functions qkifj) determine the power spectrum via 
V^{k) = hijc" /2n'^) l^fc^P, with the modes evaluated a few e-folds after fjk- Using the 
asymptotic properties of Hankel functions, the background equations and the relation be- 
tween V^{k) and 'P^{k), one arrives at well-known expression^*^. 



{k) 



2e{f]k) \ 271 J 



(5.5) 



where, on the right side, the fc-dependence is encoded in the time fik at which H and e are 
evaluated. Since H and e change very slowly during slow roll, V^lk) depends very weakly 
on k; it is almost scale-invariant. The weak k dependence is parameterized by the scalar 
spectral index ns{k), 



n 



BD 



{k) - 1 := 



d\nVf^{k) 
dhik 



-2{2e{f„) + 5{f,u)) . 



(5.6) 



(Note that in Eqs. (5.3), (5.5) and (5.6), one ignores terms which are quadratic and higher 
order in slow roll parameters.) The WMAP observations [72] provide 



Vn{K) = (2.430 ± 0.091) x 10"^ and ns{k^) = 0.968 ± 0.012 



(5.7) 



where /c* is the WMAP reference mode introduced before. For the potentials V{(j)) commonly 
employed in infiation (quadratic, quartic, exponential, etc.), the parameters e and 6 are not 
independent. For the case of a quadratic potential considered here, 6 = 0, and observational 
data (5.7) gives the values reported in (3.2) 

H{fi{k^)) = 7.83 X 10"^ and e{fi{k^)) = 8 x 10"^ 

To conclude this sub-section, let us discuss the modifications that occur if the quantum 
state of the scalar perturbation is \Q) ^ |BD). In view of the initial conditions discussed in 
section HI, 1^2) can be taken to be a 'vacuum' which is determined by a 4th adiabatic order 
basis qkiv) which differs from qf^{fi): 



qk{v) = a,qr{v) + Pkqr'{v)- 



(51 



^'^ In the cosmological literature it is common to set h ~ 1 but retain and express k in terms of the 'reduced' 
Planck mass Mp^ = H/k. 
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Here, the Bogoliubov coefficients and are time independent and are functions only of 
k := \k\. Since the basis qk{v) is also normalized, we have — \f3k\^ = 1. Physically, 
represents the number density of the BD 'excitations' with momentum k, per unit 
co-moving volume in position as well as momentum space, contained in the state The 
power spectrum and the spectral index determined by \Q) are given by 

n^{k) - 1 = n,{kr^ - 1 + . (5.9) 

If the coefficients at and are not trivial, i.e. f5k differ significantly from zero, our 
procedure will fail to be self-consistent because the background quantities extracted from 
observation will differ from (3.2), the values assumed at the start of our numerical evolution. 
We then have to seek a new solution. In the detailed numerical simulations we have carried 
out so far — with 0b > 1-15 — we find that even when the state \VL) at the onset of the slow 
roll differs from the BD vacuum, the difference is sufficiently small for the parametrization 
of the data using the BD vacuum to be adequate. 

Remark: If 0b is significantly less than 1.15, the BD parametrization will not be adequate. 
In these cases, one can use a simple 'cyclic method'. As before, one can start with the inflaton 
mass (3.1) obtained from the BD vacuum and calculate the quantum state \VL) at the onset 
of inflation. If this state differs significantly from the BD vacuum, i.e. if the Pk coefficients 
are large, one can recalculate the inflaton mass using background quantities at time ?7(/c*) 
obtained from (5.9). One can then use this value of the inflaton mass and recalculate the 
state at the onset of inflation. If this state \Vti) resulting from this ffist iteration agrees with 
\VL) we have a self consistent solution. If not, one has to continue the iteration procedure 
until there is self consistency. Of course, there is no a priori guarantee that this iterative 
procedure will converge. However, details of the pre-inflationary LQC dynamics summarized 
in section IV A suggest that not only will it converge, but the convergence may be reached 
just after the ffist iteration. Recall that: i) we are interested in only those modes which 
are observable in the CMB, and, ii) only the modes whose physical wave lengths exceed the 
curvature radius during evolution have BD excitations at the onset of inflation. As we saw in 
section IV A, this circumstance occurs only during a short interval close to the bounce (see 
the right panel in Fig. 1). At that time, the background is dominated by quantum-geometry 
effects, which are (approximately) universal, i.e., largely insensitive to the value of the mass 
m in the potential V{(j)). Therefore, after a change in the value of m as a consequence of 
the ffist iteration, one does not expect a significant change to the number of created quanta. 
Hence, the coefficients ak and Pk should remain (approximately) unchanged. This argument 
is supported by a few preliminary numerical computations for 0.93 < 0b < 1-15 where 
one obtains convergence after one iteration. These low 0b appear to exhibit a number of 
interesting features that are relevant to non-Gaussianities which will be discussed in detail 
in a separate publication. 

B. Plan of the calculations 

The scalar and tensor perturbations Q^, 7^ propagate on the quantum corrected, effective 
solution a{t) and 0(t), satisfying: 
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4 + 3i/4 + Q.^ = 0, (5.10) 

% + SH%+^% = 0, (5.11) 
where, as usual, the 'dot' denotes derivative with respect to the cosmic t, H = a/a is the 



Hubble parameter and 21 is defined in (2.14). As discussed in section III, we will assume 
that the quantum state of perturbations is a vacuum \Q) of 4th adiabatic order. We provide 
concrete examples of these states in section VC. Each of these states is determined by a 
basis which we will denote by for scalar perturbations and by for tensor perturbations. 
They satisfy (5.10) and (5.11) respectively and are normalized: 

a-^ 

euit) elit) - ekii) el{i) = % (5.12) 
Then, as discussed in section VA, the power spectra are given by 

Vrik) = h^^\e^f and Vn{k) = h^^{^^ \q^\\ (5.13) 

where the mode functions are evaluated at the end of infiation. Thus, to obtain the power 
spectra, we need to evolve the two sets of bases functions. Numerical calculations were 
carried out and we will present plots for both scalar and tensor perturbations. However, for 
brevity we will focus on the scalar modes in most of our narrative. 

To bring out the similarities and differences between the chosen 4th adiabatic order 
vacuum and the BD vacuum |BD) we will also compute the Bogoliubov coefficients Ofc 
and relating the two vacua, defined in (5.8). Using this relation, its time derivative and 
the normalization condition, it is straightforward to obtain expressions for and /S^: 

ocu = ^~a^i)[qk{i) q^^ii) - q,{i) 4^°*(t)] , 

= -t~a^i)Mi) g,^°(t) - q,ii) (t)] . (5.14) 

Note that although the right sides of (5.14) contain time-dependent terms, ak, Pk themselves 
are time independent and they satisfy the identity — |/3fcP = 1- These properties will 
be used to monitor the numerical accuracy of our simulations. 

The numerical evolutions were performed using the internal Mathematica numerical in- 
tegrators (adaptive 4*^ order Runge-Kutta). 



C. Illustrative examples of states 

As discussed in section IIIB, the quantum states \Q) of interest are vacua of 4th adiabatic 
order which satisfy an additional initial condition: At the bounce time t = 0, the renormal- 
ized energy density in the state \Q) is negligible compared to the universal energy density 
Pmax in the background. These conditions allow an infinite class of states However, in 
each numerical simulation we need to work with a specific choice. 

Our detailed simulations were carried out using the state = |0°^^) which was defined 
in section IV B. As discussed there, the state is rather simple to construct. The definition 



35 



is tied to an instant of time t and is 'local' in the sense that it is sensitive only to the 
scale factor a and its first four time derivative evaluated at t. Consider a mathematical 
example in which a is constant during small time intervals around times ti and t2- In this 
case, the prescription would lead to the standard vacuum in Minkowski space-time at these 
two times and the 'particle creation' from ti to ^2 would be exactly as in Parker's original 
work of [68, 70] which initiated the study of quantum fields in cosmological space-times. As 
a second example, consider the case in which the i derivative of = a/ a vanishes in an 
small interval around a time to- Then the state would be indistinguishable from the BD 
vacuum for high frequency modes, i.e., modes for which terms of the order 0{{k/Ha)^) can 
be neglected. Finally, as explained in section IV B, 10°*^^) has the attractive feature that, 
in a precise sense, it minimizes the uncertainties in the fundamental canonically conjugate 
fields at an instant of time. These properties motivated our use of the 'obvious vacuum', 
tailored to the bounce time, in detailed numerical calculations. 

However, to develop intuition for whether the main results are sensitive to the specific 
choice of state, we also carried out several simulations using three other states. We will 
summarize these choices both for completeness and because some of them may be useful in 
future investigations. 

• Zero energy-state at the bounce: In the absence of a larger group of isometries to 
single out a preferred quantum state, we can further constrain the family of 4th- 
order adiabatic vacua by imposing additional physical conditions. Considering the 
key role played by back-reaction in our truncation scheme, one natural requirement 
is to ask that the initial state have vanishing expectation value of the renormalized 
energy density. In the context of adiabatic regularization it is possible to construct 
4th adiabatic order states with this property at any given time [76]. These states are 
tailored to a given time, in the sense that their energy density is not zero at later 
times. Although this condition does not select a unique 4th-order vacuum, it narrows 
down considerably the possibilities. 

• Zero stress-energy state at the bounce: A stronger prescription to select a natural 
vacuum state is to demand that the state have vanishing expectation value of the 
renormalized energy-momentum tensor at a given time. When the infra-red cut-off in 
the momentum integrals is sufficiently large, this selects a unique state at the bounce 
[76]. However, for smaller cut-offs (and general times), this state does not exist. As 
in the case of the obvious 4th order adiabatic vacuum at the bounce (see subsection 
IIIB), for those cases some of the modes defining that state need to be modified. 
This introduces a freedom in the definition and the resulting state no longer has zero 
expectation value of the energy-momentum tensor. Nevertheless, this prescription is 
useful, in particular when the cut-off is sufficiently large (i.e. 0^ is sufficiently large). 

• It follows from our discussion in the beginning of section IV A that, if we re-scale tensor 
modes via Xk = a Ck, then Xk can be regarded as modes propagating in Minkowski 
space but with 'time dependent frequency' wifj). (This statement extends also to 
the scalar modes.) Following the procedure used in Minkowski space to construct 
the preferred vacuum state, one can imagine using 'instantaneous positive frequency 
modes' at a given time, with frequency w{fi). However, this naive choice needs to 
be modified both for low and high k. For low k the reason is that the 'frequency' 
w{fi) becomes imaginary. For large k the frequency is positive, however the modes do 
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not define a 4th adiabatic order state and hence, if (for example) the energy density 
is to be well defined, at large k the prescription needs to be changed. While there 
is large freedom in those modifications, if one restricts oneself to observable modes, 
the freedom becomes largely irrelevant. Therefore, this methods is well suited for 
quick calculations of power spectra and provides a useful way to analyze the effects of 
different initial conditions. 

In all these three cases, we found that the main features of the power spectra and the 
energy density in the perturbations were the same as those calculated with |0°*^^) and re- 
ported sections VD and VI. In this sense the results are robust and not tied to the choice 

\n) = |o°^^). 

D. Numerical evolution 
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FIG. 4: Power-spectrum for the scalar perturbations (crosses) for = 1.15. The solid line shows 
the average of the points and the (green) dashed line is the standard inflationary power-spectrum 
which assumes the BD vacuum at the onset of smooth roll. The deviation from the Bunch-Davies 
power-spectrum at low k is not numerical noise. Rather, modes with these frequencies experience 
background curvature during their pre- inflationary evolution and are thus excited. The state |0) 
contains excitations over the BD vacuum |BD) with quanta of these modes. This leads to a highly 
oscillatory power spectrum discussed in the text. 

Recall that the parameter space for the background is dictated by the value (pB assumed 
by the background infiaton at the bounce. Ideally one would like to explore as much of this 
parameter space as possible. However, because of the exponential relationship between 0^ 
and a(t(/c*)) (see Tab. I and Fig. 3), even a relatively modest value of (f>B results in extremely 
long integration times, which rapidly becomes computationally prohibitive. Therefore, we 
have had to restrict our attention to (pB ^ 2.0. Fortunately, as we discussed in section IV B, 
this interval covers the region of the parameter space that is physically most interesting. For, 
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our numerical simulations bear out the expectation based on those physical considerations: 
the main results are essentially insensitive to the value of 0b once 0b exceeds 1.2. 

From an observational point of view, the most significant result is the scalar power- 
spectrum, plotted in Fig. 4, for 0^ = 1.15 and the 'obvious' 4th adiabatic order vacuum, 
i.e. for \Q) = 10°*^^). We will now comment on various features of the plot. 

First, the physical argument in section IV B suggested that modes with k ^ /clqc ~ 3.21 
will not be excited because their physical wave length would be smaller than the curvature 
radius throughout the pre-infiationary evolution. Therefore, these modes would be in the 
BD vacuum at the onset of inflation. This reasoning is borne out because for large k the 
LQC power spectrum for these modes with \Q) = \0°^^) as the initial state at the bounce 
is essentially the same as that obtained in the standard inflationary scenario with |BD) as 
the state at the onset of inflation. Second, for modes with lower k values, the LQC power 
spectrum shows a highly oscillatory behavior. This phenomenon has been noted before (see 
e.g. [8, 77, 78]). Its origin can be traced back to Eq. (5.9): 

P^{k) = P|°(A;)|afc + /3fc|2 = P^^ik) {l + 2^ + 2Re {ak(3*,)) , (5.15) 

where in the second equality we have used the normalization condition — |/3fcp = 1. 
The oscillatory behavior of the power-spectra arises from the final interference term, due 
to the rapidly changing relative phase. This oscillation is so fast in k that in any realistic 
observations -which have a finite k resolution- they would be 'averaged out'. Therefore, in 
[18] we introduced 'bins' in the k space with a width 0.5ipi and averaged the oscillations in 
each bin. It turns out that the resulting plot is indistinguishable from the one obtained by 
simply neglecting the interference term in Eq. (5.15). If we do so, the LQC power spectrum 
P^{k) is given simply by rescaling P^{k) with (1 + 2|/3fcp). Since the number density of 
excitations in the mode k is given by this simplified version of Eq. (5.15) brings to 

forefront the fact that the modification of the power spectrum can be traced back directly 
to the creation of excitations by curvature. The number of excitations decay rapidly as k 
increases because modes with smaller wave lengths do not experience as much curvature 
during their pre-inflationary evolution. Hence the jik coefficients decay rapidly with k, 
damping the amplitude of the oscillations. 

The relation between the LQC and the BD power spectra is brought out more clearly in 
Fig. 5 where we plot the ratio P^{k) / P^^ (k) . To define the BD vacuum, mode functions 
were not approximated by Hankie functions for all times. Rather they were taken to be 
Hankie functions at a time when the physical frequency of the mode k^ was 10^H~^, to 
ensure that all the modes of interest were well inside the Hubble radius, and then evolved 
numerically. The circles represent the actual data points which show oscillations while the 
solid (red) curve is the plot without oscillations (or, obtained by binning as discussed above). 
Since 0^ = 1.15, it follows from Table 1 that the reference mode used in the WMAP data is 
/c* = 9.17. It is clear from Fig. 5 that, at this value, the LQC and the standard inflationary 
predictions are almost indistinguishable. Consequently, although the LQC state at the onset 
of inflation differs from the standard BD vacuum for low k, nonetheless the LQC prediction 
is in agreement with the WMAP values (5.7) of the amplitude and spectral index which are 
evaluated at the time t(/c^). 

So far we focused on the background geometry with 0b = 1.15. What is the situation with 
other values? The physical considerations of section IV B suggest that if we examine the 
excitations in any one mode k — ignoring for the moment the issue of whether it is observable 
in the CMB — then the LQC power spectrum would not be sensitive to the specific value 
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FIG. 5: Ratio of the LQC power spectrum P^{k) for scalar perturbations to the standard one, 
P^{k) which assumes the BD state at the onset of slow-roll inflation. The circles represent the 
actual ratio which exhibits rapid oscillations at low k due to a rapidly varying relative phase. The 
line is the smooth part of this ratio, which can either be produced by binning the exact data, or by 
removing the interference factor, as explained in the text. This smoothed ratio is simply 1 + 2Nk, 
where A^^^ is the number of k particles. 

of 0B- The same should hold for the ratio P^{k)/P^^{k) since the standard inflationary 
power spectrum makes no reference to (pB at all. Is this expectation borne out? Fig. 6 shows 
that the answer is in the affirmative. Specifically, we plot the fractional difference between 
the scalar power-spectra, calculated using different values of (pB- This fractional difference 
is less than ^ 0.2% over the range of k in which there is a significant deviation from the BD 
power-spectra. The true difference is, in fact, smaller than this, since the value at large k is 
dominated by numerical error (see Section VE). 

These results lead us to the following interesting overall picture. If we look at the entire 
range of wave numbers k, the ratio of the LQC and BD power spectra is essentially universal, 
i.e., insensitive to the value of 0b (at least for the range of 0b we analyzed in detail). Modes 
with k ^ /^LQC = 3.21 are in the BD vacuum at the onset of inflation while those with lower 
values of k are excited because they experience the background curvature in the Planck 
regime immediately after the bounce. What changes with 0b is the window in the k space 
that is accessible to observations in the CMB. As Table 1 shows, each 0^ determines the 
co-moving value of A;^, the reference mode used in by WMAP (consistent with our convention 
a|i-=o = The observationally relevant window is given by {ko, 2000A;o), where ko = /c*/8.58, 
and moves rapidly to the right as 0b increases. 

Next, in Fig. 7 we plot the power-spectrum of the tensor modes. Fig. 7 is very similar 
to the power spectrum of the scalar modes shown in Fig. 4, even though the evolution 
equation (5.10) for the scalar mode includes a potential 21 which is absent in the tensor 
modes. This is because from its deflnition (2.6), it is clear that the potential is proportional 
to m^, the square of the inflaton mass, and ~ 10~^^. Since the f that features in (2.6) 
ranges between and 1 and, for kinetic dominated bounces, < 10 during the entire pre- 
inflationary evolution, the effects of the external potential 21 turns out to be negligible for 
the range of fc's of interest to observations. As a consequence, the evolution of the two sets 
of bases functions is almost indistinguishable: — 1/3^"^; in the Planck regime, the same 
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FIG. 6: The ratio of the LQC power spectrum P^{k) for scalar perturbations to the standard 
one, P^{k), is insensitive to the value of (pB- Here we have plotted the fractional difference 
between ratios calculated using different values of i;^^, i-e. A(Plqc/-Pbd) = ( {PhQC / Pbt>)\^^ — 
(-PLQc/-PBD)|0g=0)( {PBD/Phq,c)\^^=^) , where cp = 1.1 and the solid (red) line is for cpB = 1-0, 
the dashed (blue) line is for (pB = 1-05 and the dotted (cyan) line is for (pB = 1-15. Note that the 
deviation from zero is almost entirely due to numerical precision, but even so the ratios are the 
same to a few parts in 10'^. 




amount of quanta are created in tensor and scalar modes. An immediate consequence is 
that the LQC ratio of the tensor to scalar power spectra, rLQc,(when averaged) is the same 
as in standard inflation 

_ 2Pj' _ 2P°P (1 + 2|gJ|') 2P°°(i-) _ 

1^ - i-i" (i + 2|/jfp) " ITW - • " " 

This result is shown in Fig. (8), where we plot the ratio rLQc/''"BD- 

To conclude, let us comment on an inflationary consistency relation. An important 
result of the standard inflationary scenario is the relation between the tensor-to-scalar ratio 
and the tensor spectral index. In slow roll inflation this relation reads tbd ~ "S'^f^! 
where the approximation indicates that terms of quadratic or higher order in the slow roll 
parameters are neglected. Since this expression holds independently of the inflaton potential 
and relates two independent observable quantities, it serves as a test of the standard scenario. 
Forthcoming observation of the effect of tensor perturbations in the CMB will provide a test 
of this relation. Does LQC modify this relation? We have seen that the tensor-to-scalar 
ratio remains unmodifled. However, the tensor spectral index is modifled because of the pre- 
inflationary evolution. The LQC tensor spectral index is obtained from the tensor power 
spectrum after averaging 

^ dlnP^ ^ dlnP^^ ^ dHl + 2\^l\') ^ ^bd _ ^ ^ rfln(l + 2|/3,^p) ^ ^^^^^^ 
Ink Ink dink * dink 



Therefore, the LQC consistency relation is given by 



, rfln(l + 2|/3PR , 
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FIG. 7: The LQC power spectrum for tensor modes. As for scalar modes, we have set (pB = 1-15 
and used the 'obvious' 4th adiabatic order vacuum at the bounce, |0) = jO"*^^) (crosses). The soUd 
hne is the average and the dashed Hne is the standard tensor power-spectrum assuming the BD at 
the onset of slow-roll inflation. 
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FIG. 8: Plot of the fractional difference between the LQC tensor to scalar ratio, tlqc and the 
usual prediction from standard inflation, tbd i-e. ('^lqc — '^bd)/'''bd • The line is the average of 
points (see discussion below Eq. (5.15)). 



41 



0.7 - 

0.6 - 

0.5 - 

0.4 - 

0.3 - 

0.2 - 



0.1 ^ ' ' ' ' ' ' ' ' ' — I 

6 8 10 12 14 16 18 20 22 24 

k 

FIG. 9: Modification of the standard inflationary consistency relation tbd ~ —8nf^, due to 
pre-inflationary dynamics. The green dashed line is —Srit, the solid black line is —8(nt — 
((iln(|afcp + |/3fcp)/(iln /c) and the red crosses are the numerically calculated tensor to scalar ratio. 

where, as before, we have averaged over the rapid oscillations. 

Fig. 9 shows that this relation is satisfied to an excellent approximation in numerical 
simulations. It shows the imprint left by the pre-infiationary dynamics which is potentially 
observable: a deviations from the standard infiationary prediction at low k's. 

E. Checks of the numerics 

The accuracy of the simulations has been checked by carrying out multiple tests. We 
verified that: 

• The norm of the evolved modes 6^(77) and qk{v) is preserved under the numerical 
evolution. 

• The Bogoluibov coefficients computed using Eq. (5.14) are time independent even 
though individual terms on the right side of this equation depend on time. 

• The relation la^p — 1/3^^ = 1 is satisfied at all times and for all k. 

• The two expressions for the power spectrum, Eq. (5.13) and Eq. (5.9) agree. 

Although at the theoretical level some of these are just identities, they are excellent tests 
for the numerics. 

As an example of convergence with respect to numerical precision, in Fig. 10 we have 
plotted the normalization condition, la^l — — 1 calculated using three increasing levels of 
numerical precision (the step size of the internal integrator being reduced until the relative 
change at that step is less than one part in 10^°, 10^^ and 10^° respectively). As can clearly 
be seen, we have already converged to better than one part in 10^ using the middle precision. 
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FIG. 10: These plots show the accuracy to which the normahzation of the modes is maintained 
throughout the numerical integration. The Bogoluibov coefficients given in Eq. (5.8) satisfy |ap — 
|/3p = 1. These plots correspond to a (cosmic) time ~ 5 x lO^ipi, for cp^ = 1.15. The left plot 
shows the normalization for precision/accuracy goals within the numerical integrator being 10^^" 
(lower jagged line) and 10~^^ (upper flat line), whilst the right plot is for precision/accuracy goals 
of 10~^^ (upper, rising line) and 10~^^ (lower flat line). Note the different scales on the two plots. 



Unless otherwise stated, all plots have been calculated with the middle precision level, with 
convergence tested using the higher precision. Note that there is a (relative) loss of precision 
at large k, which becomes important, for example, in Fig. 6 and also in the calculation of 
the renormalized energy density at late times. 

Another important check of the numerical accuracy comes from the fact that the and 
coefficients given in Eq. (5.8) are time independent. An example of the level to which our 
numerical integration maintains this is given in Fig. 11. Here the normalization condition 
is plotted (as a function of k) for various times during the evolution (with (pB = 1.15). 
A second (related) test comes from the time independence of the ratio PjI/Pbb, which is 
plotted in Fig. 12. In both cases one again sees the loss of precision for large k. However, 
for the range of k we are interested in, both quantities are conserved in time to better than 
one part in 10^. 

F. Summary 

In this section we have analyzed the effects of the pre-inflationary LQC space-time on 
the propagation of tensor and scalar perturbations. The numerical results shown in Fig. 5 
support the physical picture presented in section IV B. Pre-inflationary geometry provided 
by LQC has a significant effect the evolution of modes with low k, which have wavelengths 
of the same order or larger than the curvature scale /clqc at the bounce. However, the 
dynamics of modes with k ^ /clqc is largely insensitive to the background geometry; they 
essentially evolve as if they were in fiat space-time. 

The relevant questions is then: What is the range of co-moving k corresponding to ob- 
servable modes! As discussed in section IV B, this window depends on the value (pB oi 
the background infiaton at the bounce. The window is given by k^[^ = ko ^ A;*/8.58 and 
kma.x ~ 2000/cmin and, fc* increases with (pB (see Table 1). If 0b > 1.2 we have kmm ^ ^lqc, 
and the evolved state is indistinguishable from the BD vacuum at the onset of inflation 
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FIG. 11: Plot of the normalization condition of the Bogoluibov coefficients given in Eq. (5.8). In 
Fig. 10, we showed how this normalization converges with respect to the precision of the numerical 
integrator and here we demonstrate that this accuracy is maintained throughout the evolution. 
The various lines correspond to different time slices (from top to bottom (cosmic) time = 79, 7.9 x 
10^, 2.5 X 10^, 5.0 X 10^, 10^, 3.9 x 10^ t-p\, respectively) for = 1.15, using precision/ accuracy 
goals of 10~^^. 





FIG. 12: The ratio of the LQC power-spectrum to the BD one for scalar perturbations is a 
constant in time, since it can be expressed entirely in terms of the Bogoluibov coefficients, 
-Plqc/-Pbd (fc) = |ofc -l- As a check of our numerics, this ratio was calculated at each 

out-put time. Here we have plotted the fractional change in this ratio (relative to that cal- 
culated at (cosmic) time t ~ 1580 tpi) for a range of (cosmic) times (from bottom to top) 
t = 177, 354, 794, 3550, 6310, 12600, 25100, 50100, 10^ Planck seconds. A{Plqc/Pbd) - 
{PhQc/PBDii) - Plqc / PBD{i)) {Pb D / Plqc (t)) . The growing deviations at large k is an issue to 
do with numerical precision (see also Fig. 10 which shows the normalization of the modes). 
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for observable modes. Therefore, for 0^ ^ 1-2 the present analysis provides a justification 
for the assumption of the Bunch-Davies vacuum at the onset of inflation for both tensor 
and scalar modes (at least for the observable range of k). It is important to emphasize 
that this conclusion is non-trivial because for those values of (ps observable modes acquire 
trans-Planckian frequencies near the bounce and the standard treatment of cosmological per- 
turbation in classical space-time is not applicable. The introduction of a quantum space-time 
was essential to incorporate those modes in a controlled fashion and to explicitly describe 
their evolution starting from the deep Planck regime. 

Are there then any deviations from the standard predictions at all? The answer is in the 
affirmative for the small window in the parameter space given by (pB ^ 1-2. In this case, 
the state at the onset of inflation differs significantly from the BD vacuum, and the results 
for the spectrum of scalar and tensor perturbations differ from the standard predictions. 
For values of (ps not too far from 1.2, the modes for which the power spectrum deviates 
from standard predictions have wave numbers smaller than k-^. For these modes the WMAP 
observational error bars are large. For instance, for (pB = 1-15 we have kmin ~ 1-07 and 
modes for which deviation from standard predictions appear correspond to £ < 30 in the 
WMAP angular decomposition. As a result, the LQC predictions are compatible with 
current data for the power spectrum. However, the fact that the state differs from the 
BD vacuum for those modes leads to a deviation from the standard 'consistency relation', 
tbd ~ — 8nf°, for low k that may be seen in future observations (see Fig. 9). In addition, 
as pointed out in [10-12], the deviation from the BD vacuum at the onset of inflation can 
serve as a source of non-Gaussianity generated during inflation. Now, the presence of scalar 
BD quanta at the onset of inflation has been shown to lead to concrete effects in the CMB 
[15, 17], and in the distribution of galaxies [15, 16] which could be observed in the near 
future. If they are seen, these observations will signiflcantly constrain the initial state for 
scalar perturbations at the onset of inflation. One would be able to make speciflc predictions 
for non-Gaussianities originating in LQC, thereby opening a novel avenue to observationally 
probe the Planck regime. Because the window, (pB ^ 1-2, in the parameter space for which 
new predictions arise is narrow, should these effects be observed, they would essentially pin 
down initial conditions for the background at the bounce, making more directed analysis 
and more detailed prediction feasible. 

All these predictions were made under the assumption that our initial 'truncation ap- 
proximation' is valid; i.e., that the back reaction of the perturbations on the (quantum) 
background is negligible. The goal of the next section is to analyze if this is the case. 

VI. SELF-CONSISTENCY 

We began with the standard truncation of general relativity (coupled to a scalar fleld) 
that is used in the inflationary scenario and developed a quantum theory based on this 
truncation. A key question is whether this theory admits solutions which are consistent 
with the truncation approximation, i.e. in which the back reaction of perturbations can be 
neglected all the way from the bounce to the onset of the slow roll. Of course our conditions 
at the bounce on permissible states ip of perturbations are such that the approximation 
is satisfled initially. However, there is no a priori guarantee that it would continue to be 
satisfled during the entire evolution especially since it covers the 11 orders of magnitude in 
density and curvature. 

In the flrst part of this section, we will report on the results obtained by numerical 
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FIG. 13: Plot of the ratio of the energy density in the perturbations to that of the background for 
three values of 0b- Here the (red) + plot refers to 0b = 1-21; the (green) x to 0b = 1-22; and 
the (blue) ★ to 0b = 1-23 ). In order for the truncation scheme to be valid, we need to ensure 
pPcrt//3BG ^ 1) which is the case for 0b > 1.23. This ratio remains small not only close to the 
bounce, but throughout the evolution, as shown in Fig. 14. 



simulations using a specific quantum state of perturbations, = |0°^^). In the second 
part, we will use analytic considerations to argue that, if there is a state for which 
the back reaction is negligible, then there are infinitely many states 'close to' ip^ with this 
property. 

Our focus will be on scalar perturbations; the situation with tensor perturbations is 
completely analogous. So far, our numerical methods allow us only to put upper bounds on 
the energy density in perturbations. But fortunately this is already sufficient to establish 
the existence of a large class of consistent solutions. 



A. Numerical analysis 

Conceptual aspects of the back-reaction issue were discussed in section VI. D of [19]. It 
was shown that a necessary and sufficient condition for the back reaction of perturbations 
on the background \l/o to be negligible is that the expectation value {ipl p{x,i) \ip)rcn should 
remain negligible compared to the energy density pbg the way from the bounce to the 
onset of slow roll. 

Thus our task is to numerically evaluate the expression (2.29) of the renormalized energy 
density 



(0°'^-|p^^-^^(x,t)|0°^-)ren 



(2^ 



, ,2 (2t + A;2 
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(6.1) 
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FIG. 14: Time dependence of the ratio of the energy density in the perturbations to that of the 
background (for (pB = 1-23) over the entire evolution from the bounce (at = 0) to the onset of 
slow-roll inflation. The noisy peak at (cosmic) time ~ 5 x 10^ Planck seconds is due to numerical 
issues and in this region one should consider these points as an upper bound on the ratio. 



where, the subtraction-term is given by (2.30): 

C(«(fc,i) = *+^^^5^^ 8*3 ^ 

Recall that the infrared cut-off /cmin is given by kmin = K/S.58. Since the value of k^, is 
determined by 0b, the co- moving infrared cut-off is simply a fixed number for any given 
state \i/o of the background geometry. What would happen if one were to choose a lower 
krain as the infra-red cut-off? In this case, the results discussed below will change only in the 
value 0B above which we have self-consistent solutions of the truncated theory. Furthermore, 
the exponential relation between 0b and k^ shown in Fig. 3, implies that the required change 
in (pB will only be logarithmic in the change in /Cmin- 

Since kmin increases with 0b, one might expect that an increase in 0^ will result in a 
decrease in the expectation value of the energy density at any given particular time, simply 
because of the decrease in the domain of integration. Our numerics bear this out. The results 
are plotted in Fig. 13 for early times near the bounce and in Fig. 14 for all times from the 
bounce to the onset of slow roll. These results show, in particular, that for 0^ > 1.23 we 
have, 

<*!M=L(a < 0.05 Vf. (6.3) 

Pbg [t) 

Thus the test field approximation holds if 0b > 1.23. It may hold also for lower values of 0b 
because, as indicated in Figs. 14 and 15 and discussed below, what we control numerically 
is an upper bound on the energy density in perturbations. 
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FIG. 15: An illustration of the importance of numerical precision for high k. The plot shows the 
integrand p{k) of the energy density in perturbations for <j)B = 1.15 and at time t w 5 x lO^tpi. 
The corresponding background energy density /Obg is only ~ 10"^! The upper (red) solid line was 
evaluated with the lowest precision/accuracy goal IQ-^*^, whilst the lower (blue) line corresponds 
to precision accuracy of IQ-^^. The inset (where the y-axis is magnified by another factor of 10^) 
shows that there are delicate cancelations between positive and negative values of p{k). In the inset 
the upper (green) plot corresponds to precision/ accuracy goal of IQ-^'' while the lower (blue) to 

10-25. 

Numerically these calculations are challenging because the integrand rapidly builds up a 
fast (but decaying) oscillation in k (similar to Fig. 4). In addition, for large k, even a very 
small loss of precision can result in a significant error, since the renormalization procedure 
subtracts two diverging terms. This problem is exacerbated by the fact that there are 
times when even extremely small numerical errors (of the order of one part in 10^^) qqt^ 
comparable to the background energy density; see Fig. 15. But, typically, this is not an issue 
for the early phase where quantum gravity dominates and the background energy density 
is high, or the very late time evolution where the energy density in the perturbations is 
suppressed by the large scale factor. The problem is serious only for the intermediate times, 
t ~ 50, OOOtpi, when the calculations of the expectation value of the energy density are 
dominated by numerical noise. However we have used sufficient precision to ensure that this 
noise is always less than half a percent of the background energy density. This allows us to 
ensure that Eq. (6.3) is satisfied. But because of this, our numerical estimates should be 
considered as upper bounds. 
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B. Analytical considerations 

In section VIA we used numerics to provide an explicit example of a state for which 
our underlying truncation scheme is self-consistent. However, one would also like to know 
whether this particular example is special or if there are 'many' other states in the Hilbert 
space T-Li for which our truncation is also self-consistent. Of course, such questions have 
to be phrased with care because "Hi is infinite dimensional. In this section we will provide 
an analytical argument to establish the following result: If a state in "Hi is such that 
{ip\p\4')ren IS Very small compared to the energy density pbg in the background, then the 
same is true for every state in an open neighborhood of i/j in "Hi. Together with the 
numerical example lip) = |0°''^) given in the previous subsection, this result establishes that 
there are 'many' states in "Hi for which the truncation approximation is valid. Since our 
considerations are all analytic, in this subsection we will work with the conformal time fj of 
the dressed effective metric cjab (rather than the cosmic time t). 

Let lip) then be a 4th adiabatic order vacuum in "Hi such that. 



iv) < Pbg iv) W r] e (r/mi, r/finai) 



(6.4) 



where (r^ini, ''/final) is the dynamical range we are interested in. Now consider a second 4th 
adiabatic order vacuum \ip) in "Hi. Let qk and qk be any 4th order bases that define \ip) and 
and let the Bogoluibov transformation relating them be given by 



Qk = Oikqk + Pk ql ■ 



(6.5) 



One can easily check that the difference of the expectation values of the renormalized energy 
density with respect to and \ip) is given. 
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where we have used the relation la^p — = 1. From the previous equation, by taking 
into account that —\z\ < 'Re{z) < \z\ for any complex number we can obtain the following 
bound for Ap 
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where B± = \/3k\'^{l ± y 1/1/3^^ + 1). Finally, since {i/j) and are both 4 adiabatic 
order states, the f3k coefficients decay sufficiently fast with k so that the integrals in these 
bounds are convergent. Therefore, the bounds are finite. Furthermore, all time dependent 



quantities are continuous in t] whence the difference Ap 



(77) has an upper bound in 



the closed time interval under consideration. Since this bounds scales with Bf, , we can make 



supAp 



(77) arbitrarily small simply by rescaling the (3k by a sufficiently small constant. 



Thus we have an (open) neighborhood of the state parameterized (for each mode) by 
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in which the change in the expectation value of the energy density in the perturbations 
is also very small compared to the background for all times between the bounce and the 
onset of the slow roll. 

To summarize, in this section we have shown that if the state \l/o of the background 
geometry is chosen so that 0b > 1.23, then |\E'o) ® |0°^^) provides a self-consistent solution 
to our truncated quantum theory. Furthermore, there are infinitely many states 'close to' 
jgobv^ for which \^o) ® shares this property. Each of these states provides an admissible 
extension of the standard inflationary scenario to the Planck regime. Finally, so far our 
numerical techniques have provided only upper bounds for the renormalized energy density. 
It is conceivable that future work will show that the self consistent solutions exist also for 
lower values of 0b. 

VII. SUMMARY AND DISCUSSION 

In the standard inflationary scenario one works with quantum fields representing first 
order perturbations on a FLRW solution to Einstein's equations. One assumes that these 
quantum fields are in the BD vacuum at the onset of slow roll, works out their evolution and 
calculates the power spectra and other correlation functions at the end of inflation. The fact 
that one starts the calculations at the onset of slow roll is often justified on the basis that 
pre-inflationary dynamics will not affect any of the observable predictions. This belief stems 
from the argument that modes seen in the CMB cannot be excited during the pre-inflationary 
dynamics because, when evolved back in time starting from the onset of the slow roll, their 
physical wave lengths Aphy continue to remain within the Hubble radius £Hh all the way to 
the big bang. However, this argument is flawed on two reasons. First, what matters to the 
dynamics of these modes is the curvature radius $Kcurv = -\/6/i? determined by the Ricci 
scalar i?, and not 9^h, and the two scales are equal only during slow roll. Thus we should 
compare Aphy with 0lcurv in the pre-inflationary epochs. The second and more important 
point is that the pre-inflationary evolution should not be computed using general relativity, 
as is done in the argument given above. One has to use an appropriate quantum gravity 
theory since the two evolutions can well be very different in the Planck epoch. Then modes 
that are seen in the CMB could well have Aphy > ^cun in the pre-inflationary phase. If this 
happens, these modes would he excited and the quantum state at the onset of the slow roll 
could be quite different from the BD vacuum. Indeed, the difference could well be so large 
that the resulting power spectrum is incompatible with the amplitude and the spectral index 
observed by WMAP. In this case, that particular quantum gravity scenario would be ruled 
out. On the other hand, the differences could be more subtle: the new power spectrum for 
scalar modes could be the same but there may be departures from the standard predictions 
that involve tensor modes or higher order correlation functions of scalar modes, changing 
the conclusions on non-Gaussianities. In this case, the quantum gravity theory would have 
interesting predictions for future observational missions. Thus, as we emphasized in section 
IV, pre-inflationary dynamics can provide an avenue to confront quantum gravity theories 
with observations. 

In this paper we investigated these possibilities in the context of loop quantum gravity. 
In the standard inflationary scenario, one starts by truncating general relativity, keeping 
only the FLRW solutions (coupled to a scalar field) and first order perturbations thereon. 
Therefore, in our analysis we began in section II with the truncated phase space FTrun = 
FoXFi where Fq is the FLRW phase space of homogeneous, isotropic fields and Fi is the phase 
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space of gauge invariant first order perturbations. We then passed to the quantum theory of 
the entire FTrun- Consequently in the resulting framework tensor and scalar perturbations 
T and Q propagate on a quantum FLRW geometry encapsulated in a wave function \l/o. 
This quantum geometry is non-singular; the big bang of general relativity is replaced by 
a quantum bounce. Therefore we can directly face the trans-Planckian issues. We found 
that the existence of trans-Planckian modes is not a problem in itself. The key question 
is whether these modes contribute so much to the energy density in the perturbations that 
the back reaction of perturbations on \l/o cannot be neglected. This could easily happen 
especially in the Planck regime. If it did, our 'truncation' strategy would fail to be self- 
consistent. Therefore the challenge is to develop techniques to calculate the energy density 
in perturbations all the way from the bounce to the onset of the slow roll and compare it 
with background energy density. Finally, as we saw in section II B 2, LQC dynamics reveals 
several features in the Planck regime that are very surprising from the perspective of general 
relativity which has shaped intuition in the standard inflationary scenario. 

The underlying FLRW quantum geometry streamlines the analysis by providing the nec- 
essary control on calculations in the deep Planck regime. However, we are now faced with 
the challenge of investigating the dynamics of quantum fields T, Q on quantum geometries. 
At first this problem seems formidable. But fortunately there is a key simplification [19, 48]: 
Mathematically this evolution is completely equivalent to that of T, Q propagating on a 
dressed, effective metric Qab, which incorporates the appropriate quantum properties of the 
state This occurs because the dynamics of test quantum fields does not experience 

all the details of the probability amplitude for various FLRW metrics encapsulated in \l/o; 
it is sensitive only to a few moments of this distribution. The phenomenon is analogous to 
the propagation of light in a medium where all the complicated interactions of the Maxwell 
field with the atoms in the medium can be captured just in a few parameters such as the re- 
fractive index. This unforeseen outcome enabled us to 'lift' the adiabatic regularization and 
renormalization techniques, normally used in QFT in classical FLRW space-times, to QFT 
on quantum FLRW geometries, thereby providing the necessary control on the dynamics of 
T, Q on the quantum geometry \l/o. 

In section III we specified a class of initial conditions at the bounce for the quantum 
state \I^o ® of the combined system, \l/o describing the background and describing per- 
turbations. The permissible states \l/o are required to be sharply peaked on certain effective, 
bouncing trajectories in the phase space, which are well-understood. It turns out that, ex- 
cept for a tiny portion, (p^ < 0.93, of the allowed range of values for the background infiaton 
at the bounce, all effective trajectories encounter the slow roll compatible with the 7 year 
WMAP data sometime in their future evolution [6]. This was an unanticipated result. How- 
ever, it does not mean that the background states \l/o are generic. In fact, they are very 
special because they are required to be sharply peaked. For the states ip of perturbations, 
we imposed three conditions: i) Regularity: They be of 4th adiabatic order so that the 
expectation value {ipl p\ip)ren of the renormalized energy density operator is well defined; 
ii) Symmetry: They be invariant under spatial translations and rotations, i.e., under sym- 



Whilc the standard inflationary dynamics of tensor modes T is the same as that of a massless scalar field 
on a FLRW metric gab, that of the Mukhanov-Sasaki field Q representing scalar perturbations involves 
also an 'external potential' 21. For both modes, gab is replaced by a quantum corrected metric gab hi LQC. 
For the scalar mode, in addition, 2t is replaced by a quantum corrected 21. 
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metries of the background geometry; and, iii) Initial compatibility with truncation: At the 
bounce, {ip \ p\4')ren be negligible compared to the background energy density pbg- There is 
a large class of such states. We discussed in detail a specific example \ip) = |0°^^), called 
the 'obvious vacuum of the 4th adiabatic order', in section III, and mentioned a few other 
examples section VC. 

In sections IV and V we used numerical methods to evolve these initial conditions and 
obtained the power spectra for scalar and tensor perturbations. The key free parameter 
in this analysis is the value (pB of the background infiaton at the LQC bounce because it 
determines the effective trajectory on which \E'o is peaked. For ip, we first carried out a 
few simulations using various choices discussed in section VC to ensure robustness of the 
final results, and then focused on the choice = |0°^^) because it has several attractive 
properties. 

As explained in sections V and VI, several of the numerical simulations require very high 
precision (see, e.g.. Fig. 15). Furthermore, it turns out that the integration time increases 
very rapidly with (p^- Therefore, we had to restrict our simulations to c^b ^ 2. Fortunately, 
as summarized below, the interval 0b i$ 2 covers the physically most interesting range. 
Our numerical simulations brought out an unforeseen feature: The LQC power spectra are 
essentially insensitive to the precise value of 0^- However, in the CMB we have access only 
to a finite range of wave numbers. The 7 year WMAP data in particular covers a range 
~ {ko, 2000ko) where ko is the co-moving wave number of the mode whose physical wave 
length equals the radius £Hls of the observable universe at the surface of last scattering. In 
LQC, the value of kg is dictated by 0b (in our scale factor convention ob = 1). Therefore, in 
our analysis the window of modes (in the k space) that can be seen in the CMB is dictated 
by the value of 0b, and moves rapidly to the right along the k axis as 0b increases. 

Are there values of 0b then, for which LQC predicts deviations from the standard in- 
flationary scenario for modes in the observable window? As mentioned in the beginning of 
this section, on rather general grounds we expect that a mode would be excited in the pre- 
infiationary dynamics if Aphy > ^Hcuv at any time. A detailed analysis of the pre-infiationary 
dynamics of the LQC background geometry shows that such modes can be seen in the CMB 
only if (pB < 1-2 (see section IV B). In this case, for observable modes, the state has non- 
trivial excitations over the BD vacuum at the onset of inflation. If 0b = 1.15 for example, 
the reference mode used in WMAP corresponds to k^, = 9.17 in our conventions and as Fig. 5 
shows, the LQC predictions for the power spectrum and the spectral index at k = 9.17 are 
indistinguishable from those of standard inflation. Thus, the prediction is observationally 
viable. However, for k < 9.17 — which correspond to £ < 30 in the angular decomposition 
used by WMAP — there are significant differences between LQC and the standard inflation 
because of the excitations over the BD vacuum. Because the observational error bars for 
low i are large, both predictions are observationally viable. Furthermore, since the LQC 
excitations over the BD vacuum are the same for scalar and tensor modes, the ratio r of the 
tensor to scalar power spectra is unchanged from that in standard inflation. However, the 
standard 'consistency relation' tbd = — 8nf^ is modified. Future measurements of tensor 
modes should be able to distinguish between the two. This departure from the BD vacuum 
also has implications for the CMB and galaxy distribution [11-14] and observational tests 
for such effects have already been proposed [15-17]. Thus, there are differences between the 
LQC and the standard inflationary predictions if 0b < 1.2. 

What if 0B > 1-2? Then we would have Aphy ^ ^Hcuv throughout the pre-inflationary 
evolution for all the modes that can be seen in the CMB. Therefore, the physical reasoning 
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strongly suggests that they would not 'experience' significant curvature and would be in 
the BD vacuum at the onset of inflation. In this case, pre-inflationary dynamics would 
not modify standard inflationary predictions for observable modes. Numerical simulations 
made this semi-quantitative considerations precise in that our plots provide sharp figures for 
excitations over the BD vacuum for various values of 0b. In particular, if 0b > 1.4, we find 
that the excitations are so small that the departures predicted by LQC would be negligible 
for any of the forthcoming observation missions. This also implies that the fact that we had 
to limit our numerical simulations to 0b ^ 2 is not physically restrictive. 

To summarize, by analyzing the pre-inflationary dynamics in detail we arrived at two 
main conclusions. First, there do exist natural initial conditions at the bounce which lead to 
a completion of the standard inflationary scenario to include the quantum gravity regime. 
In this completed theory, one has a consistent evolution all the way from the deep Planck 
regime that accounts for the inhomogeneities seen in the CMB. Since the origin of the large 
scale structure can be traced back to these inhomogeneities, now one can systematically 
trace back the seeds of this structure to the quantum fluctuations of the initial state at the 
LQC bounce itself. Second, there is a narrow window in the 0b parameter space for which 
the state at the onset of inflation would not be the BD vacuum. While the LQC and 
the standard inflation predictions are both compatible with current observations, future 
observations should be able to distinguish between the two. Thus, there is a potential to 
extend the reach of observational cosmology all the way to the Planck scale. Of course, since 
the window is narrow, the 'a priori' probability of its being realized in Nature is small. 
This is compensated by the fact that, if observations are compatible with 0b being in this 
window, the initial conditions would be narrowed down tremendously, making very detailed 
calculations and predictions feasible. 

Finally, in section VI we investigated the issue of whether our truncation strategy is 
self-consistent, i.e., whether in the solutions 'X' we analyzed the energy density in 
the perturbations can be neglected compared to that in the background all the way from 
the bounce to the onset of inflation. This is a difficult issue on two accounts and to our 
knowledge it had not been analyzed in the Planck regime in any approach. The first 
difficulty is conceptual: one needs a systematic framework to compute the renormalized 
energy density. As discussed above, we were able to construct this framework using the 
well-developed adiabatic renormalization theory on cosmological space-times because of the 
exact mathematical equivalence between QFT on the LQC quantum geometries \E'o and 
QFT on dressed effective metrics gab determined by "^o- The second set of difficulties comes 
from numerics: because of rapid oscillations of the integrand of p\ip)ren in the k space, 
and because the background energy density itself diminishes rapidly as one evolves to the 
future of the bounce, one requires very high accuracy and numerical precision. We were 
not able to calculate \ p |'0)ren over the entire evolution from the bounce to the onset of 
inflation to the desired precision. However, we were able to provide an upper bound on this 
quantity. For 0^ > 1.23, these upper bounds suffice to guarantee that the back reaction of 
perturbations can indeed be ignored for the states \E^o ® -0 we considered. This suffices to 
show that there is a satisfactory extension of the standard inflationary scenario in which 
it is consistent to ignore the back reaction also in the pre-inflationary epoch. However, 
for lower values of 0b — particularly those for which the state at the onset of inflation is 
significantly different from the BD vacuum — we have not been able to demonstrate that 
the truncation scheme is self consistent. We hope to return to a detailed analysis of this 
issue in the near future. 
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We will conclude with a few remarks. 

• While this paper was focused on inflation, the underlying framework developed in [19] 
is much more general. Together, the two investigations offer some general lessons that could 
be useful also in other paradigms for the early universe. For example, we saw that, in the 
investigation of whether a mode would be dynamically excited, what matters is not the 
Hubble radius which often dominates this discussion but the curvature radius IHcurv deter- 
mined by the scalar curvature of the background. As Fig. 1 illustrates, the evolution of 
9icurv ill the early universe can be quite involved and therefore non-intuitive, especially in 
the regime in which Einstein's equations receive significant corrections. Therefore, it is im- 
portant to calculate this quantity throughout the period of interest. A second issue involves 
back reaction [79]. Since it is negligible during slow roll, the issue has often been ignored 
in the investigations of pre-infiationary dynamics as well as in discussions of alternatives to 
inflation. However, the issue is important especially in the Planck era and our numerical 
simulations showed that it is quite subtle. In particular, there is a correlation between the 
fact that modes with Aphy > ^cun are excited during evolution and the observation from 
numerics that these modes make significant contribution to the energy density in pertur- 
bations. Now, if there is sufficient inflation, then the modes that are relevant for CMB 
observations tend to have very small Aphy in the Planck regime. On the other hand, if there 
is no inflation at all, it is not easy to achieve sufficient expansion between the bounce and the 
CMB epoch for these modes to have Aphy <^ 5^curv in the Planck era. These general physical 
considerations suggest that in scenarios without inflation there is a danger that the back 
reaction may not be negligible especially in the Planck regime. Therefore it is important to 
check the consistency of flrst order perturbation theory in those cases. 

• In recent years, there have been a number of interesting investigations of possible LQC 
corrections to inflationary dynamics (for a brief summary, see section II of [19].) Some 
of these have focused on pre-inflationary dynamics studied in this paper (see especially 
[47, 80]). The distinguishing features of our analysis are: i) It is based on the general 
truncation strategy which has been successfully employed in LQG in other problems; ii) 
The issue of the initial state has been investigated in a stream-lined fashion; iii) Numerical 
simulations are better controlled and considerably more extensive. Furthermore, the physics 
behind the main flndings is understood and discussed in detail; iv) Issues of regularization 
of composite operators and renormalization of energy density have been addressed for the 
first time. They enabled us to systematically check if the back reaction can be ignored. 

• The truncation strategy used in this series of papers is motivated by two considerations: 
i) We begin with the phase space FTrun that underlies almost all investigations of the early 
universe; and, ii) The same truncation philosophy underlies other successful applications 
of LQG such as the calculation of the graviton propagator. Nonetheless, it would clearly 
be much better if we could start with a full quantum gravity theory and then descend 
to this truncation systematically. Could the final results be significantly different as far 
as observational issues are concerned? We believe that the answer is in the negative 
at least in cases where the truncation procedure can be shown to admit self-consistent 
solutions. For, in general relativity one routinely expects first order perturbations whose 
back reaction is negligible to provide an excellent approximations to the phenomenological 
predictions of the exact theory, and there is no obvious reason why the situation would be 
different in quantum gravity. As a simple example to illustrate our view of the effectiveness 
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of the truncated theory, consider the Dirac solution of the hydrogen atom. Because it 
assumes spherical symmetry prior to quantization, this truncation excludes photons from 
the beginning. Therefore, at a conceptual level, it is very incomplete. Yet, as far as 
experiments are concerned, it provides excellent approximations to answers provided by 
full quantum electrodynamics till one comes to the Lamb shift. We expect the situation 
to be similar for our truncated theory: Conceptually it is surely quite incomplete vis a vis 
full LQG, but we expect the full theory to provide only small corrections to observable effects. 

• However, even apart from the issue of full LQG, the current framework can be improved 
in a number of ways in the near future. We will complete our remarks by providing a few 
examples. 

-k Extensions: First, it would be desirable to extend the framework in several directions. 
Inclusion of a positive cosmological constant would be straightforward. While it would 
add a few conceptual twists along the lines of [31], the effect on pre-inflationary dynamics 
investigated here would be totally negligible because the phenomenological value of the 
cosmological constant is so small. Second, it would be worthwhile to extend the calculations 
to include multi-infiatons, especially since multi-infiatons have already been considered in 
LQG [81]. Are there any significant changes in the results on power spectra? A third 
extension would be to investigate cases where 0b is much larger than 2. We gave simple 
physical arguments that for 0b > 1-4 the LQG results would be indistinguishable from 
those of standard inflation as far as foreseeable observations are concerned. But it is not 
impossible that these considerations neglected to take into account some subtle feature of 
the background dynamics for large 0b. Only explicit calculations can decide. However, 
since numerical simulations become prohibitively expensive for large 0b, a combination of 
analytical approximations and numerics will be needed. Finally, our preliminary results 
indicate that, if one were to extend the detailed numerical calculations all the way to 0b = 
0.93, the lowest value that is compatible with the 7 year WMAP data according to [6], one 
would find that the state at the onset of inflation is so different from the BD vacuum that 
one has to significantly revise the value of the inflaton mass used in this paper. We have 
made a first pass at redoing the calculation using a cyclic procedure described at the end of 
section VA and plan to carry out a detailed calculation. However, a more systematic and 
efficient procedure extracted from the standard routines used in the WMAP data analysis 
would be highly desirable. 

•k Quantum theory of perturbations: To make direct contact with calculations in the 
standard inflationary scenario, we used a Fock representation for quantum states ip of per- 
turbations T, Q. This strategy is internally consistent and yields a well-defined theory. 
However, from a conceptual standpoint, it is highly desirable to use a 'polymer represen- 
tation' suggested by the LQG techniques. In cases when the energy density in the per- 
turbations is negligible compared to that in the background, we expect predictions of the 
polymer representation to reduce to those obtained here using the Fock representation. But 
if perturbations are described using the polymer representation, quantization of the entire 
truncated theory would be more firmly rooted in LQG, adding considerable conceptual co- 
herence. In addition, such an extension would also be valuable from mathematical physics 
considerations. The correspondence between QFT on quantum geometries and those 
on the dressed effective metrics Qab will continue to hold. However, LQG techniques are 
likely to suggest alternate regularization and renormalization schemes that are tailored to 
the 'polymer representation'. The adiabatic scheme used here has been carefully developed 
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within QFT on cosmological space-times over more than four decades [68, 69, 82-84]. And 
it handles the uhraviolet divergences satisfactorily. However, as we saw in section II C 2, it 
does not remove infrared divergences in the case of massless fields. Our current treatment of 
infrared issues is meant to be a physically motivated but tentative strategy. There should a 
better procedure that simultaneously handles the ultraviolet and the infrared regime. It is 
quite possible that such a scheme naturally descends from techniques that are well-adapted 
to LQG [50]. So far the literature on the relation between LQG and regularization and 
renormalization in QFT in curved space-times has remained rather general and investiga- 
tions that focus just on cosmological space-times may well provide richer and more detailed 
results. 

•k Initial conditions: In our choice of initial conditions for the background state 
in section III, we were motivated by physical considerations. The new ingredient is 
the role played by the repulsive force of origin in quantum geometry that causes the 
quantum bounce. It is significant on a scale of about lOfpi and effectively dilutes away 
the inhomogeneities. Since regions of this size at the bounce expand out to become the 
observable universe at the surface of last scattering, there is a new mechanism to achieve 
the extraordinary homogeneity that is needed in the initial state to explain the current 
large scale observations. However, so far the argument is only qualitative. It would be 
extremely helpful if this could be developed further via concrete calculations, or, evidence 
for or against this possibility can be gleaned from suitable numerical simulations. Another 
open issue is that of further narrowing down the initial conditions, especially for ip. The 
three conditions we imposed are well motivated but still allow a very large class of ip. We 
used the four choices discussed in section VC to verify that the power spectra are robust 
and then carried out the detailed numerical simulations using = 10°*^^) because there 
are several reasons in favor of this choice. Are there perhaps additional criteria that will 
single out this -or another state- essentially uniquely? It would be very helpful even to 
significantly narrow down the choices by using general physical principles. 

Finally, we would like to re-emphasize that in this series of papers we have focused only 
on the quantum gravity issues. Problems related to particle physics phenomenology and 
foundational issues related to the possible quantum to classical transition remain. 
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Appendix A: Inclusion of an inflationary potential in the underlying framework 

Since Ref. [19] was addressed to the quantum gravity audience, it focused on math- 
ematical and conceptual issues rather than phenomenological. Therefore, for simplicity 
of presentation we assumed that the scalar field is massless. In this appendix we sum- 
marize the modifications that are necessary to incorporate an inflaton potential, such as 
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V{(j)) = {l/2)m?(jp used in the main body of this paper. 

The discussion of the homogeneous phase space Fq is the same as that in [19] except that 
the scalar (or the Hamihonian) constraint So now contains the potential V{(j)) (see Eq (2.2)). 
The conceptual steps that led to the reduced phase space V of first order perturbations are 
also the same as in [19] and the treatment of tensor perturbations in [19] does not require 
any modifications. However, for scalar perturbations, expressions of the Hamiltonians of 
the final gauge invariant Mukhanov- Sasaki variables and curvature perturbations TZ^: are 
now more complicated. Our discussion of the required changes will largely follow Langlois' 
Hamiltonian treatment [85] . However, because of our primary motivation, we aim at finding 
expressions of the Hamiltonian that are well suited to loop quantization. Note also that 
the procedure uses only the constraint equations; we do not assume that the background 
space-time satisfies full Einstein's equations. 

As discussed in [19], the first order constraints generate the gauge transformations for the 
linear perturbations. One can solve these constraints and isolate gauge invariant variables 
which suffice to label the points of the reduced phase space Fi. As discussed in [19], the 
Mukhanov- Sasaki variable provides a convenient choice. To arrive at them, note first 

that, in the notation of [19], the first order metric perturbations ({abik) can be expanded 
into scalar, vector modes and tensor modes 

+ V2VP ki^aVb) + ^ Tf^ {XaXb - VaVb) + V2T^'^ {x^aVb)) (Al) 

where A; is a unit vector in the k direction and k,x,y constitutes a field of orthonormal 
triads with respect to the comoving fiducial 3-metric Qab in the momentum space. The 
Mukhanov- Sasaki variable is then given by 

where (pj: is the first order perturbation in the infiaton. 

As noted in the main text, dynamics of gauge invariant variables is governed by the part 
§2 of the second order constraint obtained by (ignoring the terms that are linear in the 
second order perturbation but) retaining terms that are quadratic in the first order fields. 
Since this Hamiltonian is gauge invariant, it can be expressed in terms of the Mukhanov- 
Sasaki variable and its conjugate momentum pl®^ 



^ I f d?k 

S2 [^'^hom = a\ = 



2 J (2vr) 



a2 



(Q) ^ 



2(^^)1%! 



(A3) 



with 



^ = -9 .4.fi_2 + 1271^—^ - — — + — — - + ^ . (A4) 



aH^n^^ aH^ n^a) d0 2 d02 



Here, we have used the lapse function tailored to the conformal time and the scale factor 
a and its conjugate momentum 7r(a) normally used in the cosmology literature (translation 
to the (i^, b) variables used in LQC can be readily carried out using Eq. (2.1)). Modulo 
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standard factor ordering ambiguities can be promoted to a well-defined operator on the 
background Hilbert space "Hq. 

Another, particularly useful form of Q arises from using §o = to eliminate the explicit 
dependence on 7r(a), 



a2 I \Ek + Ej \ Ek + E. d(h 2 d02 3 ' ' ^ ' 



where we have defined, 



E, = ^ , and E^ = a'fV (0) , (A6) 



which are closely related to the kinetic and potential energy density. Further defining the 
ratio, 

E, + E^ + aH^V {<P) 



one finds, 



= y (^f^ (0) - 2A/fK, + y,, + (^) j - y (21 + A:^) • (A8) 

It is this form of the scalar Hamiltonian that we use in Eq. (2.5). 

Finally, as explained in section VA, comparison with observations is simplest in terms 
of co-moving curvature perturbations 7?.^. To arrive at these gauge invariant variables, one 
can perform a (background dependent) canonical transformation within T^^\ generated by 

S = "-Pfn^, + Im , (A9) 

where / and z are arbitrary function of the background phase-space variables to begin with, 
but determined below. From Eq. (A9) one finds, 

7^,= ^Q,, and pf ) = ^pf + /7^, . (AlO) 

The choice of / and z is geared to simplify the form of the Hamiltonian in terms of TZj: and 
px . Therefore we will now sketch the procedure and conclude with the explicit form of z 

k 

that is needed in the expressions (AlO) of these quantities. The Hamiltonian in terms of TZf^ 
and p^ is not used in the main text but may be useful for future investigations. 

Recalling that Eq. (A9) is a 'time' dependent canonical transformation, with the evolution 
given by §0) "we find. 
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One can check that the equations of motion for TZ are independent of the choice of the 
background function /. However, we can choose it to simphfy the expression of S! 
Specifically, for 



/(7^) 



52^'' is diagonal: 
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(A12) 



[iVi 



hom 




2Q' 
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Until this point the form of the function z has not been fixed. A particular simplification 
occurs for the specific choice 



-3 P(^) 



7 Pid,) 



(A14) 



(This is the relation used in section V). Using this expression, one can explicitly evaluate 
the various Poisson brackets, for example 
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Using this in Eq. (A13) gives 
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(A16) 



This form is particularly useful for cosmological calculations, however it is clear from 
Eq. (A14) that the canonical transformation between and TZ^ is ill defined at points in 
the trajectory where p(<^) = or 7r(a) = 0. In particular, p(<^) can vanish during the evolution 
between the bounce and the onset of slow-roll inflation. Therefore, as discussed in VA, 
during this phase of evolution is well defined, whilst TZj: is not. 

Finally, as we noted at the end of section II B 1, the presence of the m^0^ potential also 
gives rise to some mathematical subtleties in the discussion of the quantum theory of the 
homogeneous background. These are summarized in section 4.3 of [7] and will be discussed 
in detail in [56]. 
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